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We present the application of the inverse scattering method to the design of semiconductor het- 
erostructures having a preset dependence of the (conduction) electrons' reflectance on the energy. 
The electron dynamics are described by either the effective mass Schrodinger, or by the (variable 
mass) BenDaniel and Duke equations. The problem of phase (re) construction for the complex 
transmission and reflection coefficients is solved by a combination of Pade approximant techniques, 
obtaining reference solutions with simple analytic properties. Reflectance-preserving transforma- 
tions allow bound state and reflection resonance management. The inverse scattering problem for 
the Schrodinger equation is solved using an algebraic approach due to Sabatier. This solution can 
be mapped unitarily onto a family of BenDaniel and Duke type equations. The boundary value 
problem for the nonlinear equation which determines the mapping is discussed in some detail. The 
chemical concentration profile of heterostructures whose self consistent potential yields the desired 
reflectance is solved completely in the case of Schrodinger dynamics and approximately for Ben- 
Daniel and Duke dynamics. The Appendix contains a brief digest of results from scattering and 
inverse scattering theory for the one-dimensional Schrodinger equation which are used in the paper. 
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I. INTRODUCTION 

A semiconductor heterostructure can be modelled by a 
system of equations describing (with a certain degree of 
completeness and precision) the state of the system. The 
equations depend on a set of structural and compositional 
data (SCD). Essentially, these are the spatial dependence 
of the chemical composition (including dopant profiles), 
the applied external fields, etc. The system's behavior 
(response) is described by functional data (FD), such as 
the electric or thermal conductance, the energy depen- 
dence of the electron transmittance, the wavelength de- 
pendence of the optical absorption coefficient, etc. The 
FD can be computed using the solution of the equations, 
and are thus functionals of the SCD. 

To design a heterostructure for a certain application is 
to find a set of SCD, which is physically (and technologi- 
cally) achievable, such that the values of a chosen subset 
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of FD will be within desirable ranges. The designer solves 
thus an inverse problem: inverting the functional depen- 
dence of the FD on the SCD. This problem is rather 
ill-posed. The desired ranges of FP may be unachiev- 
able. Generally speaking, even if a certain desired set of 
FD values is achievable, the set of SCD which achieves it 
is not unique. This absence of uniqueness is not bad in 
itself. If several solutions can be obtained, then one may 
further optimize the design in terms of other properties 
which were not included in the original specifications. 
The difficulty is mathematical. Inverting a one-to-one 
functional dependence can be a formidable task, which is 
further aggravated by the lack of uniqueness. 

A brute-force approach to the problem is always pos- 
sible: computing the FD for a set of achievable SCD 
(ideally all) and selecting the best. In practice, brute 
force optimization is restricted to rather small sets of pa- 
rameters describing the FD. This is mostly due the fact 
that any conceivable penalty function will be non-convex. 
Its graph will have a rather complex multi-valley shape. 
The search algorithm is forced to do a thorough investi- 
gation of this landscape and will eventually fail through 
run time limitations. The intelligent designer will par- 
tially avoid such restrictions by noting trends, trying to 
break the design into combining manageable blocks and 
locally improving promising configurations. 

Thus, cases when the solution of the inverse problem is 
realizable by methods which arc less costly than the brute 
force approach can be rather useful. Even if one has to 
simplify somewhat the physical model, precious insights 
on new promising configurations can be obtained. 

In this paper we will try to review the possible applica- 
tions of inverse scattering techniques to some aspects of 
heterostructure design. Our physical model for the elec- 
tron states in the heterostructure will be the (one band) 
effective mass approximation for the envelope function of 
the (conduction band) electrons. 

The simplest approach is to assume a constant effec- 
tive mass. Then, the envelopes of, electron states satisfy 
an effective Schrodinger equationu. Choosing the 2-axis 
along the growth direction, 

iftj^*(r, t) = - A- Att(r, t) + U(z)9(r, t). (1.1) 

Here, m e is the (conduction band) electron's effective 
mass. The potential is 



U(z) = £(z) + U ext (z) + $ sc (z), 



(1.2) 



where £(z) — £ C ond\c{z)\ is the (conduction) band offset 
(assumed to depend only on the local chemical composi- 
tion, c(z)); U ex t{z) is the (possibly equal to zero) external 
applied potential and we lumped in $ sc (z) the potential 
of the ionized (donor) impurities (dopants) and terms 
which will make the full potential U(z) self-consistent. 
Various models can be considered for $ sc (z): a Hartree 
self consistent potential, exchange-correlation corrections 
can be incorporated. The only requirement is that <J> sc (z) 



has to be an explicitly defined functional of the full po- 
tential U(z). 

The next step is to take into account the spatial de- 
pendence of the effective mass. The Schrodinger equa- 
tion with con stant effective mass (SE) for the envelope 
function, (IT), is replaced by the BenDaniel and DukcEl 
equation (BDD): 



2m(z) 



V*(r,t) + U{z)^(r,t). (1.3) 



Here the effective mass of the conduction band electrons 
m{z) = m con d[c{z)] is assumed to depend only on the 
local chemical compositi on c (z). The self-consistent po- 
tential is given again by ( |l.2[ ). The & sc (z) term is now a 
functional of both m(z) and U{z). 

In the following we w ill con side r only stationary states 
of the equations (1.1) and (1.3). Furthermore, since 
the potential depends only on the coordinate along the 
growth direction, z, the motion in the perpendicular 
plane is free. Setting 



(1.4) 



-Et/h) 



in the Schrodinger ( |l . l[ ) and BenDaniel and Duke fll.q ) 
equations, where E is the energy, r±_ = (x,y,0) and 
q_i = ((Zi,<Z2;0) are, respectively, the coordinates and 
components of the quasi-momcntum in the directions 
perpendicular to the growth axis, we obtain the one- 
dimensional Schrodinger, 



ip"{z) + [fc 2 



V(z)] i>{z) = 0, 



(1.5) 



and BDD 



mooip^z) 




m(z) 


+ 



m(z) 



V(z) 



ip(z)=0, (1.6) 



equations for ip(z). Here, and in the following, we use the 
notation for the electron effective mass in the em- 
bedding material: TOqq = m e for SE and = m(±oo) 
for BDD. We introduce the notations 



k = 



y/2m 00 E 



V(z) = 



2m 00 U(z) 
H 2 ' 



(1.7) 



Throughout this paper, the square root function is de- 
fined with non negative imaginary part: Im (y^j > 0. 

For real positive E in (1.7), k > 0. The prime will often 
be used for derivatives. 

The inverse spectKil, theory for the one-dimensional 
Schrodinger equatiorMa has been successfully applied to 
some optimization problems. for the bound states in semi- 
conductor quantum wellsHEI. Inverse scattering the oryllEl 
for the one-dimensional Schrodinger equation, (1.5) with 
q^ = 0, shows how one can recover the potential in ( |l.5| ) 
from the knowledge of the scattering data: the complex 
transmission and reflexion to the right/left coefficients 
{T(fc), R±(k)} for all real values of the wave-number k. 
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Widely used in electric circuit modelling00, it has been 
recently applied for designing heterostructure Bloch wave 
filter il. 

We want to solve the following problem: let the elec- 
tron dynamics b e giv en by either the SE, (1.5), or by the 
BDD equation, (|l.6| ), with the self-consistent potential 
(1.2). Find chemical composition and dopant profiles, 
going to constant limits at infinity, such that the het- 
erostructure defined by these data has a given energy 
dependence of the electron reflectance at q ± = and a 
given operating temperature: 



K(E) = \R±(k)\' 



(1.8) 



The zero of the energy scale is chosen at the conduction 
band minimum for the asymptotic composition at infin- 
ity. 

In the SE case, the mini-bands are parabolic in q ± 
and (1.8) will hold for q ± ^ with E changed to 
E + h z q 2 L /2m e . T his is no longer true in the case of 
the BDD equation (1.6). For sufficiently small \q±\, the 
mini-bands will be approximately parabolic only as long 

.*3 



as the 



m(z) 



term in the effective potential in (L6) can be 



treated as a first order perturbation . O ne could also se- 
lect a nonzero value of q 2 at which (1.8) is valid, such as 



the one corresponding to the transverse thermal energy 
at the desired operating temperature. 

Since only the energy dependence of the reflectance is 
given, the first step in solving our problem is to find the 
sets of scattering data (SD) which are compatible with 
lZ(E),i.e. find the phases of the scattering data. In sec- 
tion HI we show how to construct SD which correspond 
to real valued potentials V(z) with exponential decay at 
infinity from 1Z(E). Physically, one might expect that 
the potential is determined by the its bound state ener- 
gies and its resonances. The reflectance TZ(E) embodies 
only information on the transmission resonances: sharp 
minima of the reflectance. There is another type of res- 
onance, the reflection resonances, analogous to the res- 
onances that occur in three-dimensional potential scat- 
tering on a spherically symmetric potential. These are 
sharp phase variations of the reflection coefficients. The 
information on reflection resonances and on the bound 
states is not apparent in the reflectance. 

The transformations of the SD, which do not change 
the reflectance, will be discussed. Using these transfor- 
mations, we will define reference solutions for the recov- 
ery of the SD from the reflectance. The reference solu- 
tions have no bound states and simpler analytic proper- 
ties. Combinations of reflectance-preserving transforma- 
tions can then be used to obtain the SD of other solutions 
to the phase reconstruction problem from the reference 
solutions, by dressing them up with bound states and 
reflection resonances. ■— ■ 

We will use Pade approximation methodalj to repre- 
sent the scattering data and to find parameterizations 
for a large class of solutions, corresponding to poten- 
tials which tend exponentially to zero at infinity. We will 



find that on this type of input data, the (re) construction 
process amounts essentially to finding the roots of some 
polynomials and grouping them into subsets. 



In section 



[II we present a simple and efficient algo- 



rithm for solving the inverse scattering problem for scat- 
tering data in the form obtained in sectiordp.— We will 
use results due to Kay, Moses and SabatierEirtl for the 
inverse problem with rational coefficientsEEI. 

In section IV we present the variable mass,, unitary 
mapping of the SE to the BDD equatiorM^lil We 
formulate the boundary value problem which must be 
solved for determining the coordinate transformation 
which defines the mapping, given the material relation 
between the effective mass and the band offset. This 
ill-conditioned problem can be solved by a a shooting 
method. In the case when the relation between the mass 
and the offset is linear, the solution takes a simpler form. 
We also give an efficient perturbative method for solving 
the mapping equation. 

In section [v| we show that obtaining chemical compo- 
sition profiles and self-consistent potentials for them in 
the SE in the inverse scattering approach is simpler than 
obtaining the self-consistent potential for a given chem- 
ical composition profile. We also discuss the functional 
equation which must be solved for obtaining the chemical 
(effective mass) profile corresponding to self-consistent 
potentials in the BenDaniel and Duke's equation. 

For the reader's convenience, in Appendix ^ we give a 
brief outline of results from scattering and inverse scat- 
tering theory, which are needed and often referred to in 
the main body of the paper. 



II. PHASE RECONSTRUCTION 

Inverse scattering theo ry for the one-dimensional 
Schrodinger equation, (|l.5[) with q ± = 0, 



i/>"(z) + [k 2 - V(z)] i>(z) = 0, 



(2.1) 



on which we give a primer in Appendix |a|, shows that 
one can recov er a fast decaying and piecewise continuous 
potential in (2.1) from the knowledge of the scattering 
data (SD): the complex transmission and reflexion to the 
right/left coefficients {T(k), R±{k)} for all real values of 
the wave-number k if there are no bound states. If bound 
states are present, knowledge of the SD is not sufficient 
for unique recovery of the potential. If the number of 
bound states is exactly n, then a n-parameter family of 
potentials gives exactly the same scattering data. 

The SD are completely determined by one of the re- 
flexion coefficients and the values of the energies of the 
bound states. The phase of T(k) can be obtained from 
a logarithmic dispersion relation, (see e.g. the book by 
Chadan and Sabatienl3, XVII.1.5 ). T he other reflection 
coefficient can be obtained from ( All ). 

In this section we assume that the reflectance, TZ(E), is 
known on the positive energy half-axis, E > 0. We want 
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to construct sets of SD which satisfy ( |l.S| ) . Since the val- 
ues of the SD for scattering by short-range and piecswisej 
continuous potentials must satisfy the constraintstL3'Eil 
which are enumerated at the end of Appendix A 2 , the 
function TZ(E) cannot be arbitrary. It must be non- 
negative and smaller than— unity, with the exception of 
TZ(0), which is genericallycHl equal to 1. If the potential 
is piecewise continuous, TZ{E) must go to zero no slower 
than E~ 2 for large values of E. 

The transmittance T{E) = |T(fc)| 2 , where T(fc) is 
the complex transmission coefficient, is readily recovered 
from 1Z(E) +T(E) = 1. Thus, we know the absolute 
values of the scattering data and we need the phases. 



A. Reflectance-preserving transformations 

The problem of finding the phases of the scattering co- 
efficients knowing only their absolute values on the real 
axis is underdetermined and has an infinite number of 
solutions. Before considering the phase (re)construction 
problem, we will introduce two types of transformations 
which modify the phases of the scattering data without 
changing the reflectance. 

Let the set of scattering data 



{T(k), R+(k) R-(k)}, 



(2.2) 



be a solution of the phase recon stru ction problem, i.e. 
the scattering coefficients satisfy (|l.8|) an d the conditions 
enumerated at the end of A ppendix [A3. Then, as men- 
tioned above, if the SD (2.2), has no bound states, there 
is an unique solutio n to the inverse scattering problem: 
a potential V(z) in (2.1) such that the SD calculated for 
this equation coincide with ( |2.2| ). If T(k) has n simple 
imaginary poles in the upper half plane, that is n bound 
states, then a n-parameter family of potentials can be 
constructed, such that the SD of each potential coincides 
with (U). 

Let now A > be a positive number, such that k = iA 
is not a pole of T (k), i.e. that E = — h 2 \ 2 /2m e is not 
a bound state of (2.2). Define a new set of SD by the 
transformation: 

T(fc) 
R±(k) 



(fc + iA) 
(k - iA) 
(iA + fc) 
(iA - fc) 



T(fc), 



R±(k), 



(2.3) 
(2.4) 



The new set of SD will also satisfy the conditions set out 
at the end of Appendix A 2, so that a (n + l)-parameter 
family of potentials can be constructed with each poten- 
tial having the SD (OJ27 



2.4) can also remove bound 
([2.2|), have a bound state for 



The transformation (f 
states. If the initial SD, 

E = —H 2 \ 2 /2m e , then the transformation (2.3-2.4) with 
A changed into —A in the right-hand sides of the equa- 
tions, transforms that bound state into an anti-bound 



state - an imaginary pole of T(k) in Im(fc) < 0, leaving 
all the others in place. 

Let C, be an arbitrary complex number with nonzero 
real and imaginary parts. Then, we can define a second 
t ype of reflectance preserving transformation of the SD 

& 



T(fc) 
R-(k) 

R+{k) 



T(fc), 

(fc-C)(fc + C) 
(k-C*)(k + Q 

(fc-C*)(fc + C) 
(fc-C)(fc + C*) 



JZ_(fc), 



R+(k). 



(2.5) 
(2.6) 

(2.7) 



Here and in the following, we use the nota tion * fo r com- 
plex conjugation. The transformed SD, (^^2jj), have 
the same reflectance and bound states as ( [2.2 ) and sat- 
isfy the conditions enumerated at the end of Appendix 
A 2 . Using the inverse scattering method, one can con- 



struct from the SD, (2.5-2.7), a new n-parameter family 
of piecewise continuous potentials which goes to zero at 
infinity. 



The transformation (2.5-2.7) has a simple interpre- 
tation. Assume that the imaginary part of £ is much 
smaller than its real part, |Im(£)| <C |i?e(£)|, and the ini- 
tial reflection coefficients, R±(k), are slowly varying on 
the scale |Im(£)| near fc = ±Re(C). Then, the new scat- 
tering data ( |2 . 5) - |2 . 7[ ) have a reflection resonance of width 
|Im(C)| at fc = ±Re(£). Indeed, the phases of the new 
reflexion coefficients vary by ±2tt in a small interval of 
width 2|Im(C)| centered on fc = ±Re(£). 

A third type of transformation adds purely imaginary 
resonances to the reflection coefficients leaving the trans- 
mission coefficient unchanged: 



T(k) 
R+(k) 

R-(k) 



T(fc), 
(iA + fc) 



(iA 
(iA 



fc) 
k) 



(iA + fc) 



R+(k), 



R_(k), 



(2.8) 
(2.9) 

(2-10) 



where A is an arbitrary real number. The new SD will 
also satisfy the conditions enumerated at the end of Ap- 
pendix A 2. The reconstruction of the potential is done 



exactly as for the second type of transformation. 

A sequence of transformations of the second and third 
type can be written as: 



T(fc) 
i?_(fc) 



T(fc), 
S N (k) 
S N (-k) 



R-(k), 



(2.11) 
(2-12) 

(2-13) 



where S^ik) is any polynomial whose zeros are invariant 
with respect to reflection through the imaginary axis, 
i.e. if £ is a zero, then — £* is also a zero of the poly- 
nomial. Such polynomials, normalized by the condition 
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SW(O) = 1, satisfy 5jv(fe) = [Sn(— k)]* for real k. Let us 
recall that an arbitrary polynomial of degree n, n„(a;), 
with n„(0) = 1, can be expressed through its zeros: 



n 



X 

x, 



(2.14) 



Here Xi, i = 1, . . . ,n are the zeros of Tl n (x) (including 
multiple ones according to their algebraic multiplicity). 
Thus, SW(fc) is completely determined. 



B. The reference solutions 

Let us assume that we have found a solution of the 
phase reconstruction problem. Gcnerically, it will have 
some bound states and reflection res onances . Using suit- 
ably chosen transformations of type (2.3-[2.4|) one can ob- 
tain from it a solution for which the transmission coeffi- 
cient has no poles in the upper complex half-plane (i.e. 
no bound sta tes) . Th en, by a sequence of transforma- 
tions of type (fO-2.10) with suitably chosen parameters 
one can find a solution for which the reflection to the left 
coefficient has no poles or zeros in the upper half plane. 

We will call this solution of the phase reconstruction 
problem the left reference solution. The left reference 
solution's transmission coefficient, T( r ~)(fc), and the re- 
flection to the left coefficient, R^L \k), are analytic and 
have no zeros in Im(fc) > 0. In a similar way, we can de- 
fine the right reference solution, for which T^ r+ \k) and 



R^ +> (k), are analytic and have no zeros in Im(fc) > 0. 
The left and right refer ence soluti ons are connected by a 
transformation of type ( |2.1l| - pTl~3 ): 



T {r+ \k) =T (r -\k) =T«(fc), 



R { : + \k) 



~>N+M 



(A) 



(-*) 



R ( r\k), 



R 



(r+)/,x _ S N+M (-k) ( r _) 



(k) 



SN+M{k) 



R^'ik). 



Here, the polynomial 



S N+M (k) = A N (k)B M (-k), 



(2.15) 
(2.16) 

(2.17) 



(2.18) 



where the zeros of the polynomial Aj^(k) / B]ij{k) co- 
incide (including multiplicities) with the poles/zeros of 



R^ \k) in Im(fc) > 0. Taking into account (pTj) the 



polynomials are completely determined. 

The reference solutions are in a certain sense the max- 
imally non-symmetric solutions. Indeed, as shown in 
the next section III, the potentials corresponding to the 



left /right reference solutions are identically zero for x < 
/ x > 0. Other solutions of the phase reconstruc- 
tion problem can be obtained from the reference ones 
by adding bound states and reflection resonance s wit h 
reflectance- preserving transformations of type (p^- 2.10 ). 



The left reference solution's R^ \k) is analytic and 
has no zeros in Im(fc) > 0. Thus, the logarithm 

ln[i2_ ^(fc)] is also analytic in Im(fc) > 0. The phase 
of the left reference solution (which is equal to the imag- 
inary part of ln[i?L r_) (fc)]) can be obtained from the log- 
arithm of its absolute value (which is equal to the real 

part of ln[i?^ \k)]), usins-a— (subtracted) logarithmic 
dispersion relation (See e.(?.E3'Eil). 

We will proceed in a different manner, which is more 
adequate with the physical context. 

The effective mass approximation is valid only for 
energies within an interval not exceeding several hun- 
dred milliclcctronvolts (meV) near the T point minimum 
in the Al x Ga\- x As system (or in the lattice matched 
Ini-x-yAl x Ga y As systems). Thus, two sets of scatter- 
ing data having the same (or close) low and intermediate- 
energy behavior, but whose exact high-energy behavior is 
different, can be considered equivalent. We need a good 
approximation of the SD in the physically relevant range 
of energies, which obeys the high-energy constraints set 
forward-in the Appendix |A|. The Pade approximation 
methodic is a good framework for that. An added bonus, 
which will be apparent in the following section III , is the 



simplification of the calculations needed for recovering 
the potential. 



C. Pade phase reconstruction 

We start with approximating the input design data 
for the reflectance by a type II [p,p+q+2] Pade 
approximantlij: 



K(E) = 



-Pp(E) 

Qp+q+2{E) 



(2.19) 



where V P (E) and Q p + q +2(E) are polynomials of degrees 
p and, respectively, p + q + 2 with p, q > 0. Since the 
reflectance must be non-negative and less than 1 for all 
E > 0, the polynomials must satisfy 



< V P (E) < Q 



p+q+2 



(E), 



(2.20) 



on the positive half-axis. The second inequality ( 2.20| ) 
is strict for all E > and becomes an equality only for 
E = 0. This-ensures that 7?.(0) = 1 as it should be in the 
generic caseE3. We can rewrite it in the form 



Q p+q+2 (E) =V P {E) + E1C 



P+9+1 



(E), 



(2.21) 



where the polynomial K. p + q +i(E) > for E > 0. The 
polynomial V P (E) is normalized by setting 



V P (0) = 1. 



(2.22) 



Finally, the transmittance is approximated by the 
[p+q+2, p+q+2] Pade approximant: 



T(E) = 1 - K(E) = 



EK 



p+q+l 



(E) 



V V {E) + EK p+q+1 (E) 



(2.23) 
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We can choose the coefficients of the polynomials in 
( [2. 19 ) as the parameters. These can be obtained from 
standa rd ty pe II Pade fitting routines. Taking into ac- 
count (|2.20D with equality at E = and (|2.22|) , the fit is 
obtained by solving a system of 2p+g+2 linear equations 
with 2p + q-\-2 unknowns (the polynomials' coefficients), 
which make the fit exact at 2p + q + 2 chosen points. 
Thus , the [p,p+q+2] Pade ( 2.19| ) is fully determined. Eq. 



( |2.23 ) does not introduce additional parameters. 

As we will see further on, it is advantage ous t o repa- 
ramctcrize in terms of the zeros and poles of (2.19). Since 
the reflectance and the transmitta nce a re real and non- 
negative, the zeros and poles of (2.19) and ( 2.23| ) are 
either real negative or come in complex conjugate pairs. 
The only exception to this rule are eventual real positive 
zeros of V p {E\ which have an even order of degeneracy 
(generically —2). The heterostructure is transparent to 
Bloch waves (maximal transmission resonances) at the 
energies corresponding to these degenerate zeros. 

Now, we wan t to find the left reference solution by 
solving eq. (|h8|), 



R { r-\k) =K(E), 



(2.24) 



with the reflectance given by ( 2.19|) . The reflection to the 
left coefficient, R^~\k), will be sought as a [p,p+q+2] 
Pade approximant in the variable k 



R (r- )(k) ^_P E (k) 



Q 



P+q+2 



with the normalization 



P P (0) = Q P+9+2 (0) = 1, 



(2.25) 



(2.26) 



which agrees with (2.22). As me ntione d above, \k) 



for the left reference solution of ( 2.24 ) is analytic in the 
upper half-plane and has no zeros there. Then, the zeros 



of both the denominator and numerator of ( 2.25 ) must 
lie in the lower complex half-plane. A furth er constraint 
on the roots follows from the relation ( A13 ), 



[R-{k)\* = R-(-k). 



(2.27) 



This relation holds only if for each r, which is a zero/pole 
of ( 2.25| ), — r* is also a zero/pole. 

Let us now use (2.14) t o rep resent all the polynomials 



involved in the equation (2.24) as products 



nr=i(i-w(i+w n?=i(i-£M) 



n^? +2 (i 



k/q ){l + k/ qj ) 



nS +2 d 



(2.28) 



Here pi and p^, i — 1, . . . ,p are the zeros of P p {k), re- 
spectively V P (E), while qj and q_y, j = 1, . . . , p + q + 2 
are the zeros of Q p+q+ 2(k), respectively Q p+q+ 2(E). 



Taking into account the relation (1.7), E = h 2 k 2 /m e , 
we find the following relations between the zeros of the 
polynomials: 



1j 



h 2 pf 

2m e 
H 2 q 2 
2m e 



i=l,. 



J = l,...,p + q + 2. 



(2.29) 



(2.30) 



These relations solve the problem up to the ambiguity 
of the signs of the square roots. Let us show that the 
reference solution is unique. 

As mentioned in the Introduction, the square-root 
function as maps the complex plane cut along [0, +oo) 
onto the upper complex half-plane, i.e. if E is in the cut 

plane, then Im (v^) > 0. Since the zeros and poles of 

the reference solution cannot lie in the wpper half-plane 



the sign choices in (2.29-2.30) are uniqueE 

tiPi = -\/2m ( $^ i = l, 
Hqj = - 



\/2m e (\j\ j = l, 



,P- 



q 



2. 



(2.31) 
(2.32) 



It remains to show that if r is one of the zeros, respec- 
tively poles, of (2.25), then — r* is-also a zero, respectively 
pole. This is obvious if r is realE^I or imaginary. Other- 
wise, both the real and the imaginary par ts of r 2 are 
nonzero. Then, r* 2 is also a zero (pole) of ( [2.19 ), since 
as noted above these come in complex-conjugate pairs. 
But our definition of the square-root function leads to 

(Vr*) 



Thus, ( 2.25 ) with the zeros and poles given by ( 2.31 



2.32) is the unique left reference solution for the reflec- 
tion to the left coefficie nt if the reflectance is given by 
the Pade approximant ( [2. IS ). We can readily recover 
the complex transmission coefficient of the reference so- 
lutions using the same approach. 

The (complex) transmission coefficient of the reference 
solution is given by the [p+q+2,p+q-|-2] (diagonal) Pade 
approximant 



T^(fc) 



kK t 



p+q+l 



) p+q+2 (k) 



(2.33) 



normalized by setting the coefficient of k p+q+1 in the 
polynomial K p+q+1 (k) equal to the coefficient of k p+q+2 
in the already determined polynomial Q p+q+ 2(k), so that 
limfe^oo T(k) — 1. Let and t,-, j — 1, . . . ,p + q + 1 be 
the zeros of the transmittance, (2.23). 



r P (t) = q 



p+q+2 



(«)■ 



(2.34) 



Repeating the reasoning that led to ( 2.31 
the zeros Kj of the polynomial K p+q 



we recover 



■y/2m e l 



3' 



and 



K, 



P+q+l 



(k) 



3 



n 



,p + q + l, 



P+Q+1 {Ki 



k) 



"9+2 

=1 yj 



(2.35) 



(2.36) 
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Finally, using the relation ( A16| ), we recover the left 



reference solution's \k) as a [2p+q+l,2p+2q+3] 
Pade approximant: 



R 



(k) = 



K 



p+q+l 



(k) P p (-k) 



p+q+l 



(-k) Q p+q+2 {k) 



(2.37) 



The first fraction in the above expression is a phase fac- 
tor. 



Thus , if the reflectance is given in the Pade form (2.19), 
the left reference solution is uniquely determined. The 
right/left reflection coefficients of the right r eferen ce so- 
lution are given by right-hand side of ( [2.25 )/( 2.37 ). The 
potentials corresponding to the right /left reference solu- 
tions are mirror images: V^ r ~\x) — V^ r+ \— x). 

Other solutions can be obtained by using transforma- 
tions of type (f2~3| - [Q|) to introduce bound states int o th e 



reference solutions. Then, transformations of type (2.5 



2.13) can be used for reflection resonance management. 
A distinguished class of solutions, having the same trans- 
mission coefficient as the reference ones and involving no 
additional parameters, can be obtained by redistributing 
reflection resonances between the reflection coefficients to 
the left /right. This is achieved by factorizing the poly- 
nomials P p and K p+q+ i 

P P (k) = P p+ (k)P p _(k); K p+q+1 (k) = K n+ (k)K n _(k), 

(2.38) 

into f actor s which satisfy the complex conjugation rela- 
tion ( 2.27 ) for real k. Here, p + + p_ = p, n + + n_ = 
p+q+l and all the factors are n ormalized by P p± (0) = 1. 
Applying the transformation ( |2.1l[ - pTT^ ) with S(k) — 
K n _ (k)P p+ (—k) to the left reference solution, we obtain a 
solution h aving the same transmission coefficient T^ r \k), 
given by ( |2.33| ), and the reflection coefficients are given 
by the [n±+p,n±+p+q+2] Pade approximants 



R_(k) 
R + (k) 



K n _(k) P p+ (-k)P p _(k) 
K n _(-k) Q p+q +2{k) 

K n+ (k) P p+ {k)P p _(-k) 
K n+ (-k) Q p+q+ 2{k) 



(2.39) 
(2.40) 



III. SOLUTION OF THE INVERSE SCATTERING 
PROBLEM FOR THE SCHRODINGER 
EQUATION WHEN THE SCATTERING DATA 
ARE RATIONAL FUNCTIONS 

In this section we will present the solution of 
the inverse scattering problem for the one-dimensional 
Schrodinger equation in the case when the scattering 
data are given in the Pade approximant form we ob- 
tained in the preceding section [5 C| . The inverse scatter- 
ing problem in the case of rational coefficients has been 



first considered by Kay aad-JVIose Significant re- 

sults are ,dtt6,to SabatierHil We will follow Sabatier's 
approacht§ll3 quite closely. 



As shown in the Appendix A 2, the potential in (2.1) 
can be recovered from the transformation kernels. 



V{x) = +2—K±(x;x + 0). 
ax 



(3.1) 



The transforming kernels K±(x ; y ) are the solutions of 
the Marchenko equations, ( A33 ), ( A35 ) 



K_ (x; y)+M-(x + y)= / dsM_ (y + s)K- (x; s) = 0, 

J —oc 

(3.2) 

K+{x- y) + M+(x + y)= dsM+(y + s)K+{x; s) = 0. 

J X 

(3.3) 



The Marchenko kernels, M±(u), are given by ( A34 ) 



M ± (u) = - 



+ OG 



dKe ±iKU R±(K) 



E(^r 2 < 



(3.4) 



Here, Ct are the normalization constants of the bound 



states, defined in (LA14 ). 

Inspection of ([3.2]) and (3.3) shows that the first vari- 
able, x, enters the equations only as a parameter. Also, 



for negative/positive x, the integral in (3.2) /( 3.3) in- 
volves only negative/positive values of s. Since we need 
only the contact values of the transformation kernels in 
( |3.ip , it is natural to use the K_ version of ( 3.1 ), obtained 
by solving ( ft.2j ) for negative x, and the solution of (3.3) 
for positive x. Takin g in to account the definition of the 
Marchenko kernels, ( |3.4| ), we see that we can close the 
integration contour into Im(fc) > in the expression of 
the relevant kernel (M + /M_ for x > 0/x < 0). In partic- 
ular, if one of the reflection coefficients is analytic there, 
its contribution to the corresponding Marchenko kernel 
is identically zero. In the absence of bound states, this 
implies that the corresponding M± = and, therefore, 
the potential will be zero on the corresponding half-axis. 

If the reflection coefficients are rational functions, like 
the solutions of t he p hase reconstruction problem we ob- 
tained in section II C , they have at most a finite number 
of poles in Im(fc) > 0. The corresponding Marchenko 
kernels (3.4) will be given by finite sums of exponentials, 
which go to zero at the corresponding infinity. For sim- 
plicity's sake, we will assume that all the poles are simple. 
The multiple pole case can be dealt with as a limiting 
case of pole confluence. Then, the corresponding kernel, 
M + (u)/M_(u) for u > 0/u < 0, is 



M±(t*)=i 



(3.5) 



jen± 



Here the sets of (all distinct) complex numbers \i>j } - en 
consist of the (simple) poles of R±(k) in Im(fc) > and, if 
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bound states with energies —h 2 \ 2 /2m e are present, 
The coefficients are either the residues of the corre- 
sponding reflection coefficients or, if ir^ comes from a 
bound state, — i (Cp) 2 . 



Substituting the separable Ma rchenko kernels (3.5) 
into the Marchenko equation (3.2) for x < 0, we obtain 



(3.6) 



where 



Yr( x )= e - iv i x + &yK^{x-y)e-' w i y . (3.7) 



Multiplying ( p.6[ ) by e -1 ""^ and integrating with respect 
to y we obtain a system of linear equations for 

the $:(f2_) unknowns Y~(a;). Here we used the notation 
#(J7) for the number of e lements in the set £1. Substitut- 
ing the solution into (3.6), we obtain after a little algebra 



K^(x;x-0) 



dx 



-Trln 



1 - c 



D"e 



-ixN 



(3.8) 



Here, N and D are #(il_) square matrices with 



N mj = 5 mjV j ; 



(3.9) 
(3.10) 



Solving in the same manner the equation for K + {x\ y), 
we obtain 



K + (x;x + Q) = —Trln 
dx 



1 - e 



(3.11) 



where the elements of the #(f2+) square matrices N + and 
D + (x) are 



(3.12) 
(3.13) 



We can substitute now Q and (pl| ) into (§!]) to 
obtain 



V(z) = -8Tr 



(3.14) 



1-0*0) 1-Q ± (a;) 
where the ± signs are for x > 0/x < and 

Q±(^) = e 2i l-l N± D ± . (3.15) 



Thus, the potential is given by (3.14) in terms of a 
trace of the inverse of finite matrices. 



IV. VARIABLE-MASS MAPPING: 
SCHRODINGER'S TO BENDANIEL AND 
DUKE'S EQUATION 



In section III we have presented a simple algorithm 
for solving the inverse scattering problem for the one- 
dimensional Schrodinger equation in the case of rational 
SD. 

If the conduction electron dynamics is described 
by the BDD equation, we will not start by posing 
the inverse scattering problem for that equation from 
scratch. Instead of |*hf4[— }ve w ^ use a f amn y °f um ~ 
tary transformational3'B3't3. The transformations map 
the one-dimensional Schrodinger equation, ( |1.5| ) with 
q ± = 0, and potential Vs(z), 



^"( z ) + [k 2 -Vs(z)}i>(z) = 0, 



(4.1) 



into BenDaniel and Duke equations, ([1.6]), with q ± = 0, 
and potential Vbdd(z), 



m(z) 



+ [k 2 -Vbdd{z)]^(z) = 0, (4.2) 



with a variable effective mass m(z) and a potential 
Vbdd{z). The effective mass and the new potential are 
functionally related to the parameters defining the uni- 
tary transformation. Since the mapping is unitary, the 
SD f or the equation (4.2) will be identical with the SD for 
l| ) . We start thus with a Schrodinger reference equa- 
tion, fl4.ip , which is the solution of the inverse scattering 
problem. From it we obtain a family of BDD equations 
with the same scattering data. The problem is to choose 
among these transformations those which map the solu- 
tion of the inverse problem for the SE onto acceptable 
BDD equations. 



A. Unitary mapping 

Let us introduce a (nonlinear) coordinate transforma- 
tion 



z = X(x), 



(4.3) 



which maps the interval (— oo,+oo) into (— oo,+oo). 
Here the function X(x) is a smooth monotonically in- 
creasing function: X'(x) > 0. The monotonicity ensures 
that a unique inverse transformation exists: x = Z(z), 
with 



Z[X(x)] = x 



(4.4) 



and X[Z(z)] — z. The inverse function Z(z) is also 
smooth and monotonically increasing, with 

Z\z) = l/X'[Z(z)}. (4.5) 

Let us associate with the coordinate transformation 



(4.3) a mapping Ux of the space of square-integrable 
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functions, L 2 , which transforms each element of the space 
/ 6 L 2 into U x f with 



(U x /)(aO = yOC{x)f[X(x)]. 



(4.6) 



Making the change of variables (4.3) in the normalization 
integral 



(/I/) 



+ 00 



Az\f{z)\ 2 = dxX'(x)\f[X(x)}\ 2 

J -co 

Ux/|0x/) = (/|U x Ux|/), (4.7) 



we see that the mapping Ux maps the square integrable 
functions into square integrable functions conserving the 
norm, i.e. Ux is isometric on L 2 . Since the inverse trans- 
formation. 



(V/) [z) = yfmj£)S\Z{z)\. 



(4.8) 



exists and is non-singular, the mapping Ux is unitary: 



U 



x 



U 



x 



(4.9) 



We want to see the effect of the transformatio n Ux 
on functions satisfying the Schrodinger equation (4.1). 
We will consider only the coordinate transformations for 
which the function X(x) is twice diffcrentiable with a 
piecewise cont inuo us third derivativeEEl, X'" . Substitut- 
ing (fO) into (4.1), we obtain after some algebra, 



X'(x) 



[X'(x)} 2 
Here \{ x ) — {^x^){x) and 

X"'{x) 



+ [k 2 -Wx{x)]x{x) =0. (4.10) 



W x {x) =V s [X(x)] 



5[X"(xW 

2[X'{x)} 3 4LY'(z)] 4 



(4.11) 



The equation ( 4. 1C ) satisfied by the t rans formed func- 
tion xi x ) resembles the BDD equation ( |4.2| ) if we set 



(4.12) 



m(x) = moo = m cond [c(x)]. 

The function X(x) must satisfy 

lim X'(x) = 1. 

x — ► ioo 



(4.13) 



This ensures that m(x) tends at infinity to the constant 
limit moo. We will also require that X'{x) — 1 decays at 
infinity faster than | a; | 2 5 for some 8 > 0. In this case, 
the limits of X(x) — x at ±oo are finite. This ensures that 
the Jost solutions of the original Schrodinger equation are 
mapped into Jost-type solutions of the transformed one. 

Let f±(x; k) be the Jost solutions of (4.1), which obey 
the boundary conditions (A3): 



lim c Tikx f±(x;k) = 1. 

x — *±oo 



(4.14) 



Then, applying to them the mapping ( |4.6| ), we obtain 
near the corresponding infinities 



lim e* ikx (U x f±)(x;k) 



Jkd± 



(4.15) 



Here, 



r-\-oo rO 

d+= &t [X'(0 - 1] ; d_ = / d£ [X'(0 - 1] ; 

JO J -co 

(4.16) 

where we assumed ^(0) = 0. Let f^\x; k ) be the Jost- 
type solutions of the transformed e quati on (4.10), defined 
by the same boundary conditions (4.14). Taking into ac- 
count the asymptotic beh avior s of the Jost functions of 
the Schrodinger equation (LI), we find the relation be- 
tween the scattering data for the transformed equation 
and those of the original one: 



T^ x \k) =c ikid +- d -)T(k), 



R±(k). 



R { ± ] (k) =e 



±2\kd± 



(4.17) 
(4.18) 



Thus, the reflectances of the original and the trans- 
formed equations are equal. 

In the following subsection we will study the variable 
effective mass mapping nonlinear differential equation in 
the general case of an arbitrary dependence of the band 
offset on the concentration. We will show that the physi- 
cally realizable concentration profiles for devices embed- 
ded in a material of homogeneous composition can be 
obtained by solving a (nonlinear) boundary value prob- 
lem. A necessary condition for the existence of achiev- 
able solutions is that the dependence of the conduction 
band offset £ c on the effective mass be a non- decreasing 
one. If this condition is not satisfied, it may be still pos- 
sible to embed the device in a periodic super lattice, a 
case which will not be discussed in this paper. In section 
1VD we consider in more detail the case when £ c (m) is 
a linear function, an approximation which seems quite 
reasonable in the AlGaAs system. In this case, the solu- 
tion of the nonlinear boundary value problem can be then 
expressed through the canonical solutions of the Cauchy 
(initial value) problem for a third-order linear differential 
equation. We also discuss methods for obtaining stable 
approximate solutions. 



B. Differential equation 



Let us now compare the p oten tials Wx, (4.11), and 
Vbdd- The latter is given by (1.2) and (1.7): 



Vbdd(x) 



2m a 
19 



[£(x) - £ (oo) + $ sc (z)] • (4-19) 



Here, the band offset is 

£ (x) = £cond[c{x)}, 



(4.20) 
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and we set U ex t(x) = 0. The zero of the energy scale is 
chosen at the value of the band offset for the embedding 
(asymptotic) composition. 

The band offset, £ (x) and the effective mass, m{x) de- 
pend on the position only through th eir de pende nce o n 
the local chemical composition c(x) 1 ( 4.12 ) and ( 4.20 ). 
We will consider only the case when the function m(c) is 
invertible. Then, we can substitute the inverse function 
c(m) into £ con d{c): 



£cond{c) = £ c (m). 



(4.21) 



In this section we consider only unbiased structures 
with negligi ble el ectron density. Then, in the square 
brackets in ( 4.19 ), & sc ( x) is equal to zero. Equ ating 
the two potentials ( 4.1l[ ) and (4.19) an using (4.12), we 
obtain a third-order nonlinear differential equation for 
X(x): 



X'"(x) 5[X"(x)} 2 



Vs[X(x)] 



(4.22) 



2[X'{x)f 4[X'(x)] 4 

= {£ c [m 00 X' 2 (x)] - £ c (tooo)} • 

The values of X'(x) are restricted to the physically 
achievable interval 



moo V m r 



(4.23) 



Here m and m are the minimal and, respectively, the 
maximal values of the effective mass, m(c), in the physi- 
cally achievable chemical composition range. 

The equation (4.22) does not depend explicitly on x. 
Setting 



X' = S [X(x) 



(4.24) 



substituting this into (4.22) and replacing the derivatives 
with respect to x according to 



_d_ _ dX__d_ _ d 
fa~~dx~dX~ cLY' 

we obtain a second-order equation for S, 



(4.25) 



2^ - 3^2- = ^ [£ c (mooS 2 ) - fc(rnoo) - U(X)] 



S 



S 2 



r, 2 



(4.26) 



Here S' is the derivative of S with respect to X, U(X) = 
H 2 Vs(X)/2m O0 is the pote ntial in energy units. The 
asymptotic condition (4.13) reads now S(±oo) = 1. 



C. Boundary value problem 

Let us consider potentials V(X) which are identically 
equal to zero outside some interv al \X- , X+], For X out- 
side this interval, the equation ( 4.26| ) does not depend 



explicitly on X so that its order may be further reduced. 
Setting 



and 



m(X) = rriooS 2 
S'(X) = Q(m(X)), 



the equation ( [4.26|) becomes 
^Q 2 



dm 



3^ = — 

m h 2 



£ c (m) - f c (moo) 



(4.27) 
(4.28) 

(4.29) 



The solution of ( |4.29| ) is 

4m 5 



Q z (m) = 



n 2 



£ c {y) - f c (moo)]M 2 d^. (4.30) 



If m is close to moo, the right hand side of (4.30) goes 
to zero as (m — moo) 2 . Indeed, for small |/i — m.oo|, the 
term in square brackets in (4.3C) is approximately equal 
to £' c (moo) (m - moo) and 



Q\ m) „ 2£ 'c(™°°) (m _ moQ f + [(m- moo) 3 ] 



h 2 



(4.31) 



Since Q 2 is non-negative, inspection of (4.31) shows 
that the effective mass may tend to a constant limit at 
infinity if and only if the band offset dependence on the 
effective mass, £ c (m), is a non- decreasing function. 



d£ c (m) 
dm 



> 



(4.32) 



In plain words, if the effective mass does not follow the 
band offset, then the potential V(X) cannot be embed- 
ded in an alloy of homogeneous composition. Noting that 
in this case the device may be embedded in a periodic 
superlattice, we will restrict ourselves here to the case 
(p~32[). 



Let us note that equation ( [4.26|) may be solved in 
quadratures on the in terva ls (-c o. X ) and (X + , +oo). 
Indeed, s ubst ituting ( 4.30| ) into ( 4.28 ) and taking into 
account (4.27) we get 



dS 
dX 



±sign(l - S)Q{ moo S 2 



(4.33) 



where Q(m) is the non-negative square root of the right 
hand side of ( 4.30| ) and the sign is positive on (X_|_, +oo) 
and negative on (— oo, X_). 

The solution on (X +1 +oo) is given in parametric form 

by 



X = X + + siga[l-S(X+)} f 

Js 

x = x++sign[l- S(X+)} 

Js + (x + ) o-Qim^a 1 



S+(X+) Q(mooCT 2 )' 

dcr 



(4.34) 
(4.35) 
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where 35+ is the value of x which maps to X + : X + = 
X(x+). A similar representation is valid for the interval 
(-oc,X_). 

Thus, we have obtained the solutions which are regular 
at infinity. There, S{X) = 1 + [e~ K|x| ] tends expo- 
nentially to the constant limit 1 . The asymptotic rate of 
decay k may be obtained from ( 4.31 ): 



2m^ 



V2^(Woo). 



(4.36) 



Let us now briefly discuss the asy mpto tic behavior at 
infinity for the regular solutions of (4.26) for potentials 
V(X) whic h go to zero at infinity. One may readily see 
from (4.26) that S(X) must still converge to 1 at infin- 
ity. As long as the potential V(X) decays at infinity 
faster than e _K ' x ' the properties of the solutions dis- 
cussed above remain asymptotically valid. For potentials 
with slower falloff at infinity, the regular solution tends 
asymptotically to the regular solution of the equation 
S" = 8 moc £' c ( moc )S/h 2 + V(X)/2. 

It is important to note that any regular solution of 
(4.26) satisfies at all points X± outside the support of 
the potential the bo unda ry conditions ( 4.33 ). To con- 
struct a solution of (4.26) which is regular on the whole 
axis X, we will chose a value for the embedding alloy ef- 
fective mass, moo. The regular solutions on the interval 
(—oo, X_) can be parameterized by the value S(X- ) = t. 
The value of the derivative at X- is given by Q4.33D with 
the minus sign: 



S"(X_;r) = -sign [1 - r] Q(m 00 r 2 



(4.37) 



where Q(m) is the positive square root of the right-hand 
side of (f£|g). 

Then, we integrate the equation ( 4.26 ) numerically 



with the initial values defined above over the support 
of the potential up to X + , obtaining S(X + ;t) and 
S' (X+ ; t). If the solution is regular then it must satisfy 
( 4.33 ) with the plus sign at X + : 



S'(X+;t) = sign [1 - r] Q [ moo S 2 (X+; r)] . (4.38) 



If they exist, the solutions of equation ( }4.38| ) give the 
values of r for which we c an find regular solution of the 
differential equation (4.26) on the whole real axis. To be 
achievab le th ese solutions must also satisfy the physical 
bounds ( fj~25| ). 

Numerically, the shooting method outlined here is 



rather ill-conditioned. Since the general solution of ( |4.26| ) 
is singular, multiple-precision arithmetic has to be used 
for the integration of the differential equation if the sup- 
port of the potential is not short enough. In the following 
section we will examine in more detail the important case 
when the dependence of the band offset on the effective 
mass is linear. We will show that in this case the solutions 



of the nonlinear equation (4.26) can be found among the 
solutions of a third-order linear equation. 



D. Solution of the boundary value problem in the 
case of linear dependence of the band offset on the 
effective mass 

For some alloys like Al c Ga\- c As in the concentration 
range < c < .45 the conduction band minimum is 
at the center of the Brillouin zone (r) and the depen- 
dence of m con d(c) and £ C ond(c) on the concentration c 
approximately linear. In the Al c Ga\- c As system, for 
< c < .45, the offset from the position at GaAs is 



£ c {fn) = B(m — m). 
The constant B is 

AE 

B = — = 9.41eWm , 

m — m 



(4.39) 



(4.40) 



Here m = .067mo is the conduction band effective mass 
for GaAs; to = .104m and AE are, respectively, the 
mass and band offset for Al A5 Ga, 55 As. m is the elec- 
tron mass. 

Substituting this into ( 4.26j ) we obtain 



25" S - 3(S'f = n^S 2 [S 2 - 1 - v(X)] 



Here. 



*& = SBm^/h 2 , 
v(X) = 4V(X)//& = U(X)/(Bm 00 ), 



(4.41) 



(4.42) 
(4.43) 



and U(X) = ?i 2 V(X)/2m O0 is the potential measured in 
energy units. 

Let us define a new unknown function 



T(X) = 3-3/S(X). 



(4.44) 



Substituting ( 4.44 ) into ( 4.41 ) we obtain the equation 
satisfied by T(X): 



{1- [l + v]T 2 } =0. 



(4.45) 



Here and whenever it docs not lead to ambiguities we will 
omit the arguments of the functions. Taking the deriva- 
tive of ( p~4E| ) we obtain a third-order linear equation for 
T: 



rpfi 



k 2 oo (1+v)T'-IkI o v'T = 0, 



(4.46) 



An arbitrary solution of (|4.46| ) will satisfy the second- 
order nonlinear equation ( [4.45D with the zero in the right 
hand side replaced by some constant JC: 



2TT" - T' 2 -«£,[! + v] T 2 = K. 



(4.47) 



If K = — K 2 ^, the solution satisfies also ( 4.45| ). 

Let F a j3 1 (X] X_) be the solution of the initial value 
(Cauchy) problem for the linear equation (4.46) satisfy- 
ing the initial conditions 
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F Lf3~f( X - 
F af3~/ ( X - 



X-) 
X-) 
X) 



■ a, 
■■0, 
: 1, 



(4.48) 



at the point A_ , whic h is a point of continuity of the 
potential v(X). Since ( 4.46 ) is linear, 

(4.49) 



F a (3 1 = CtFioo + 0Fqio + 7^001- 



The funct ion F a g y (X; XJ) also satisfies the second-order 
equation (1.47) with 



[!+«(*-)] 



(4.50) 



in th e righ t hand side. 

In IV C we outlined the shooting method for solving 
the boundary value probl em w hich leads to the regular 
(acceptable) solutions of (4.26). It is well known that 
shooting methods are prone to numerical instabilities 
even for linear boundary value problems. Another un- 
pleasant feature is the fact that we have to integrate the 
equation numerically over the support of the potential 
for each value of r. 

Let us state the boundary value problem for the func- 
tion T(X), assuming the potential v(X) to be identically 
zero outside the interval (X_, X + ). The solutions regular 
on (— oo,A_) and (X+,+oo) are 



T-{X) = 1 
T+{X) = 1 



[T(XJ) - 1] e K ~ {x - x -\ 
[T(X+) - 1] e - K ~( x -*+). 



(4.51) 
(4.52) 



Thus, the boundary conditions at X± which must be sat- 
isfied by the regular solution are 



T'(XJ) = k 00 [T(X_)-1], 
T'{XJ) = - Koo [T(X + ) - 1] . 



(4.53) 
(4.54) 



Thus, we can parameterize the regular solutions (4.51) 
by the value of T(XJ) = a. Then, T'(XJ) is given by 
(4.53). Then we integrate (4.45) up to X+, where T and 



its derivative must satisfy ( 1.54 ), whence the acceptable 
values of a are determined. 

Now, we can use the linear equation satisfied by T to 
express the solution of the initial value problem through 
the canonical soluti ons fi pp, fpip and -Fooi of the initial 
value problem for ( 4.46 ). Using ( 4.45 ) with v(Xj = 



to find T"(XJ), we find the solution satisfying (4.53) 



T a (X) = aF 100 (X;XJ) + Koo (a - l)F QW (X; XJ) 

+ K 2 oo (a-l)F 001 (A;A_). (4.55) 

At X+, T a {X) must satisfy the boundary condition 
(4.54). Whence we find a: 



a = 1 



F[ 00 (X+; XJ) + Koo [F 100 (X+;X.) - 1] 
F[ K k2 (X + ;XJ)+ Koo F lKooR 2jX+;XJ) : 



where F a p^ was defined in (1.48). Thus, the solution of 
t he b oundary value problem for the nonlinear equation 
(1.45) is expressed through the canonical solutions of th e 
initial value problem at A_ for the linear equation ( 1.46 ). 

The value of A_ can be chosen arbitrarily as long as 
the potential v(XJ) = 0. We may safely assume that 
v(X) = in a neighborhood of A_. The functions i*ioo 
and Fi K K 2 are also solutions of the second-order nonlin- 
ear equation (1.47) with JC determined from (1.5C). We 



have found the value of a, which determines the solution 
of the boundary value problem, through the solutions of 
two second order (albeit nonlinear) equations. Although 
the integration must be performed only once from A_ to 
X+, it is still numerically unstable. 



E. Approximate solutions 

As mentioned above the numerical integration of the 
differential equation over typical device lengths 20-40nm, 
having the order of magnitude of the electron mean free 
path, is rather ill conditioned. Taking the value s of the 
parameters for the AlGaAs system, above (1.4C), and a 
value TOqo « .lnioo for the effective mass of the embed- 
ding material, we obtain that the natural length scale for 
the differential equation (4.41) is 1/kqo ~ ■<Snm, which is 
comparable to th e lat tice period. 

The equation (4.41) can be rewritten as 



1 - v(X) 



2S"S - 3(S'f 



(4.57) 



Since the right-hand side of ( 4.57 ) has the small factor 
kJ 2 , one is tempted to proceed in a a naive "quasi- 
classical" way and neglect the right-hand side entirely. 
Then, 



S(X) « ^l + v(X), 



(4.58) 



and finding the coordinate transformation X(x) reduces 
to a simple quadrature. This gives surprisingly reason- 
able results. 

Finding corrections to ( 4.58)) seems a rather tedious 



task, especially for potentials w hich h ave discontinuities. 
We will note instead that from (1.43), 



v(X) 



U(X) _ U(X) 
Brrioo ~ .946^' 



(4.59) 



for the AlGaAs system. Since typical potential values 
are ilOO^OOmeU, the potential v(X) (measured in the 
natural units of the problem) is small compared to 1. A 
"small potential" perturbative approach to solving the 
boundary value problem is thus indicated. 

Let us assume that v ~ a and seek S as a series 

S(X) = 5 (0) (A) + aS (1) (X) + a 2 S {2) (X) + . . . , (4.60) 



(4.56) substitute v — * av and (4.41), rewritten as 
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S" = 35' 2 /2S + k^S [S 2 - 1 - v] /2 , 



(4.61) 



F. Example 



and expand into a power series in a. Then, after equating 
the terms with the same power of a and setting a = 1, 
we obtain a hierarchy of linear differential equations for 
the functions Sm (X), I = 0, 1, . . .. 

The first equation from (p60|) is S'LXX) = 0. The 



solution must go to 1 when X goes to infinity so that 
5 (0) (X) -> 1, while all the other S {e) (X) -> as 
X — > ±oo. Then, 



S (0) (X) = 1. 
For £ > 1, the hierarchy has the form 

q/I 2 c _ 2 

^(i) - ^oo^Ci) — ~2 K °° t '' 
3 1 

5" 2 ) _ K L 5 '(2) = ~5(i) + x ftoo 5(i) [ 3 5(i) 



2 

5(") - ^L5(£) = ^ (X; 5(!), . 



'(f-i) 



,v) ; 



(4.62) 

(4.63) 
(4.64) 

(4.65) 



Taking into account the boundary conditions for X — > 
±oo, the solution of the £-th equation in the hierarchy is 

1 r°° 

%W = -^— / d*e-"-"*-«l 

•^"■OO J —OG 

(£;5(i)i---,5(*-i)>«) , (4.66) 



which can be verified by direct substitution. 

The large Kqo limit of the first terms in the perturba- 
tive expansion ( 4.6C| ) is 



S {Q) (X) + S (1) (X) = 1 + ^e-^ x ~M0 

i (4.67) 



1+2^), 



which coincides with the first terms in th e per turbative 
expansion of the " quasi-classical solution" ( 1.5S ) . Several 
terms we checked also have this property. 

Having found S(X), we can now find the coordinate 
transformation by a simple quadrature: 



x 



dY 
o 5(F)' 



(4.68) 



In fact, this yields the inverse transformation, (4.5). Fi- 
nally, the mass profile in the BDD equation is given in 
parametric form by (4.68) and 



m(X)=m 00 S(X). (4.69) 
The concentration profile can now be readily recovered. 



We will illustrate the method by finding the specifica- 
tions for a filter with two narrow transmission resonances, 
centered at E\ < E^, with the high-energy resonance sig- 
nificantly narrower than the low-energy onetHI. Appar- 
ently, this is difficult to achieve by the direct methodEZl. 

We will start by constructing the Pade a ppro ximant 
for the reflectance. We will chose q = 1 in ( ^.19| ). This 
ensu res that the potential in the Schrodinger equation 
(4.1) will be continuous. W e wil l try to keep p, the de- 
gree of the denominator in (2.19), ass small as possible. 
Choosing the reflectance to be equal to zero at the en- 
ergies Ei and E2, we find that we need at least p = 4. 
This means two double real zeros at E\ and E^. Thus, 
the denominator is 

V 4 {E) = A 2 {E)B 2 {E) = (1 - £/£i) 2 (l - E/E 2 f. 

(4.70) 

There is some freedom in choosing the numerator, 
Qt{E). Baring other constraints, it is convenient to note 
that 0.7(E) has to have a pair of complex conjugate zeros 
in close vicinity to each of transmittance maxima E\ and 
E%. We will choose 

Q 7 (E) = [A 2 {E) + E / F{\[B 2 {E) + E / F 2 ]Q Z {E). (4.71) 

Here F\ (F 2 ) must be much larger than E\ (E 2 ), to en- 
sure that the zeros of the factors in the square brackets 
are close to the corresponding zeros of the denumerator. 
Indeed, in these conditions, the two poles near E\ are 
given by 



E k Ei [l±iy/E 1 /F 1 



(4.72) 



and a similar expression for the poles near E 2 . It remains 



to determine the positive polynomial Qz{E) in (4.71). 

Since we want sharp transmittance maxima near Ei 
and E 2 , we want that the background reflectance, 



TZ b (E) - 1/Q 3 (E), 



(4.73) 



will be close to 1 and relatively slowly varying in the 
range of energies 0, E, with E > E 2 . This can be achieved 
by choosing 



Qs(E) 



1+ i 







1 + 


(¥-01 







(4.74) 



Then, if 6 <C 1, the background reflectance is monoton- 
ically decreasing with TZb (0) = 1 and IZb(E) ps 1 — S. 
It has a horizontal inflection point at E — E/2 with 
K b (E) w 1 - 8/2. 

Neglecting 5 with respect to 1, the full widths at half 
maximum (FWHM) of the transmittance maxima are, 
respectively, 2 y /Ef/F 1 and 2- s /E\jF % . Thus, F x and F 2 
can be obtained from the FWHM of the corresponding 
resonance. 
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We chose E\ = AOmeV and E 2 = lOOmeV, with 
FWHMs equal to 9meV and, respectively, 3meV as ini- 
tial data. Proceeding as explained in sections O, pT| 



IV D we obtained an effective-mass profile which yielded 
the same reflectance as the input one. This continu- 
ous profile was digitized manually into 12 steps with 
lenghths which are integer numbers of lattice constants 
and heighths which are a combination of three concen- 
trations, as described in Table 1. The Al concentration 
of the embedding alloy is 11.4%. In Fig. 1 we present 
the transmittance calculated for the resulting configura- 
tion. The digitization has introduced a few artefacts (the 
shoulder of the low energy line and a low amplitude broad 
maxima, one of which is visible in Fig. 1). The maxima 
were broadened and shifted a little from the design data. 
Nevertheless, the ratio of FWHM of the high/low energy 
maxima of transmittance is better than 2:1. 



V. DETERMINATION OF CHEMICAL 
COMPOSITION AND DOPANT 
CONCENTRATION PROFILES OF 
HETEROSTRUCTURES WITH PRESET 
REFLECTANCE 



In the section [II we have shown how, given rational ex- 
pressions for the scattering data (SD), one can construct 
the piecewise continuous potential V(x), decaying expo- 
nentially at +oo. The one-dimensional Schrodinger with 
V(x) has the given SD. In section^ we mapped unitarily 
the one-dimensional Schrodinger equation onto a family 
of BenDaniel and Duke type equations depending on the 
function X(x) which defines the transformation. We also 
discussed in some detail the solution of the nonlinear dif- 
ferential equation, which determines the effective mass 
profile in the case of an undoped heterostructure, with 
vanishing density of conduction electrons. 

Now, we want to deal with the case when a nonzero 
densit y of conduction electrons is present and the poten- 
tials, (jT^) in the SE, ([□]), or in the BDD equation, Q, 
are self-consistent. Let the external potential U ex t{z) — 
so that the potential in the BDD equation is 



U (z) = £cond[c{z)] + $ sc (z; to; U), 



(5.1) 



where we evidenced the functional dependence of 
Q S c(z; m; U) on the effective mass profile m(z) and the 
full potential energy U(z). In the SE case the effective 
mass is constant and $ sc depends only on U(z). Inspec- 



tion of ( |5.l| ) suggests the idea that for a given full po- 
tential U(z), we can determine the chemical composition 
profile by mo ving the selfconsistent potential to the left 
hand side of (5.1). We will see that in the SE case this 
is relatively easy to do at a given operating temperature 
T a , The situation is trickier in the case of BDD dynamics, 
where we will present a perturbative approach to solving 
the functional equ ation which replaces the nonlinear dif- 
ferential equation ( 4.22| ). 



Let us start by considering the Hartree approxima- 
tion for the self-consistent electrostatic potential. Then, 
& sc {z; m; U) is the solution of Poisson's equation 



(5.2) 



Here — e is the electron charge, e(z) is the dielectric con- 
stant, e(z) = s[c(z)], and p c h{z) is the full charge density, 
the difference between the ionized donor charge density 
and the electron one. 



Pch(z) = e [n d (z) - n el (z)] 



(5.3) 



The density of donor dopant ions is made of a uniform 
background density rib and the local variation of the den- 
sity of donors in the heterostructure ni{z), which goes to 
zero for large \z\: 



n d (z) = n b + n e (z). 



(5.4) 



We assume that the doping and temperature are such 
that all the donors are ionized and that the density of 
holes is negligible compared to the density of donors. 
Then, we can neglect the valence bands. Otherwise, a 
multi-band treatment is needed. 

Th e electron dynamics is descr ibed either by the SE, 
(1.1), or by the BDD equation, (1.3). The equilibrium 
electron density, n e i(z), can be calculated from the den- 
sity of states for the corresponding equation, ( lO ) or 



n e i(z) 



+ OC 



dE v{z, E) 

1 + ^(E-v) ■ 



(5.5) 



where (3 is the inverse temperature (in energy units) and 
/i is the chemical potential of the electrons. 

A necessary condition for the stability of the system is 
that the full charge is equal to zero, 



dzp ch (z) = 0. 



(5.6) 



In particular, the limiting value of electron density at in- 
finity equals the background ion density, n e ;(±oo) = n&. 
The limiting values of the potential at ±oo will be equal 
only if dipolar moment of the charge density is zero, 



dzzpch(z) = 0. 



(5.7) 



Thus, z 2 p c h (z) must go to zero at infinity and the solu- 
tion of (5.2) is 

(5.8) 



/ + OC n. 
dup ch (u) / 
-oc J u 



e(v) 



For position- independ ent e , which we will consider in the 
case of SE dynamics, (5.8) becomes the well-known 



2ire f + °° 
$ sc (z; U) = — - I du\z - u\p ch (u) 



(5.9) 



We will consider separately the cases when the elec- 
tron dynamics is described by the Schrodinger eq uati on, 
(1.1), and by the BenDaniel and Duke equation, (1.3). 
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A. Schrodinger's equation 

To compute <& sc {z; U) using (j5j|), we need the electron 
density of states in (|5.5|) . Since the transverse degrees 



of freedom separate, ([1.4-1.5), we can integrate over the 
transverse quasimomenta and express the three dimen- 
sional density of states us(z,E) through the density of 
states Uo(z,E) of the one-dimensional SE (4.1): 



m e [°° 
u s{z,E) = —^\ dnv (z,E 
nnr Jo 



rj). 



(5.10) 



The one dimensional density of states Vq{z, E) is pro- 
portional to the imaginary part of the Green function 



Gq(z, z'\ E) of the one-dimensional SE (4.1 



vo{z 1 E) = limlm [Go(z, z; E + iS] 

7T <510 



(5.11) 



This allows a standardized treatment of the bound and 
continuum states. 

For all E with lm(E) ^ 0, the Green function is the 
solution of the equation 



Since the selfconsistent potential varies slowly with the 
temperature, at temperatures which not too far from T , 
the selfconsistent potential of the structure, which can be 
computed by standard means, will be close toJJ{z). An 
important result of scattering theory, see e.<7.£2l, is that 
the SD and the bound states of the SE are continuous 
functionals of the potential. Thus, the SD and bound 
states in an interval of temperatures near T Q will be close 
to those at T D . 

Now, a big advantage of the inverse method becomes 
apparent. In the direct approach we give the composition 
profile c(z) and the dopant profile ni{z), and have to com- 
pute the selfconsistent potential. This is a rather time- 
consumming iterative process, which has to be repeated 
after each modification of the chemical and dopant pro- 
files. Only after this are we able to obtain the scatter- 
ing data and check if they are desirable. In the inverse 
approach, we start from desirable SD and can obtain a 
whole class of specifications, parameterized by the local 
dopant profile ng, which yield the desired properties ex- 
actly at T a and approximately in some neighborhood of 
T . 



E 



h 2 d 2 

2m P dz 2 



-U(z) 



G (z,z';E) = S(z- z'), (5.12) 



B. BenDaniel and Duke's equation 



which is continuous at z = z' and goes to zero for 
z — > ±oo. In terms of the Jost functions, defined in 
Appendix A 1 , the solution is 



G (z,z';E) = r J^flf + ( z; k)Mz';k); z > z' , 



ikh 2 

G(z',z;E); z < z 1 . 



(5.13) 



Here k = \/2m e E/h and T(k) is the transmission coeffi- 
cient. 

Substituting ( pT0| ) and ( |5.1l| ) into (|J) we obtain the 
electron density 



riei(z) 



irfih 2 



AE In 



1 + e 



/3(/x-E) 



u (z,E), (5.14) 



Now, let us fix the operating temperature T = l//3 . 
The local dopant concentration profile, n^(z), will be cho- 
sen such that the conditions ( |5.6| ) and ( |5.7| ) hold a t T Q . 
This can be done in many ways. For each ne(z), (5.1) 
yields a band offset profile 



£ CO nd[c(z)] = U(z) - $ sc (z; U). 



If the values of the local doping profile ne(z) and the con- 
centration profile, c(z), which can be determined from 
(5.15) are physically achievable, a heterostructure with 
these specifications can be manufactured. Assuming 
Schrodinger dynamics for the conduction electrons, the 
selfconsistent potential at the operating temperature T a 
will be U(z). At T , the SD and the bound states of the 
electrons in the heterostructure will be those determined 
by U(z). 



In section |y| we have shown that to each monoton- 
ically increasing on (— oo,oo) function X(x), one can 
associate a unitary operator, (ix, which maps the SE 
onto a BDD type equation. There we considered only 
undoped heterostructures without occupied conduction 
electron states. Given the material relation between the 
band offset and the effective mass, we have shown how 
to solve the nonlinear differential equation which deter- 
mines the transformation X(x) and obtain the band off- 
set (and effective mass) profile which has the same SD as 
the initial SE. 



As in section IV B, we assume that the effective mass' 



dependence on the local chemical composition is invert- 
ible. Then, material relations e xist i n terms of the effec- 
tive mass for the band offset, (4.39), and the dielectric 
constant, 



e{z) = e[m(z)]. 



(5.16) 



The electrostatic part o f th e selfconsistent potential, 
<f> sc (z; to; U), is given by (5.8) in terms of e[m(z)] and 
the full charge density p c h(z; to; U). We can now repeat 
the cal culat ion of the electron density of states like in 
section V A . The transverse degrees of freedom separate 



again and, after performing the angular integral over the 
transverse quasi-momenta 



The electron density is given again by (5.E) where the 
density of states is the integral over the transverse mo- 
menta of the one dimensional density of states. S etting 
q 2 ^ = 2m ao ri/h 2 and performing the angular integral over 
the transverse quasi-momenta we obtain 
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Tfl f°° 

'i z i E ) = I dfjiy n {z,E- 

7™ Jo 



V), 



(5.17) 



where v v {z,E) is the densit y o f states of the one- 
dimensional BDD equation (4.2) with the potential 
Vbdd replaced by V n (z) = 2m O0 U ri (z)/h 2 where 



U v (z) = Ubdd(z) + t] 



1 - 



m(z) 



(5.18) 



The density of states for the one dimensional BDD equa- 
tion can again be expressed in terms of its Green function, 

vJz,E) = --limIm\GJz,z;E + iS)} . (5.19) 

7T <5i0 

For all E with lm(E) ^ 0, the Green function G n (z, z'\ E) 
is the solution of the equation 



E 



dz 2m(z) dz 



- U n (z) 



G rl (z,z';E) = S{z-z'), 
(5.20) 



which is continuous for z = z' and goes to zero for 
z — ► ±oo. Like in the case of the SE, the one dimen- 
sional BDD Green function can be expressed in ter ms o f 
the Jost solutio ns, f± (z'; k; rf) of the BDD equation ( |4.2| ) 
with p ote ntial ( |5.18] ) which satisfy the boundary condi- 
tions (A3). Here, k = V2to 00 _E/ ^. W e will write here 
only the contact value needed in ( |5.19] ) 



G v (z,z;E) 



T(k) 



ih 2 k 



f + (z;k; V )f_(z;k; V ). (5.21) 



Substituting ( 5.17 ) into (5.E) and making the change 
of variable E — > E + r\ in the double integral, we obtain 



n e i{z) 



dE 

x. J0 



d?/ 



Vr,(z,E) 

1 + e /3(B+r,-M) ' 



(5.22) 



In the SE case, discussed above in subsection |V A| , 
m(z) — rriao = m e and v v (z,E) does not depend on 
j]. This allowed us to perform explicitly the integration 
with res pect to the energy of transverse motion, r), and 
obtain (5.14). Now, the mini-bands are not parabolic, 
and fn{z, E) depends on r\. Nevertheless this dependence 
is weak for sufficiently small rj. Let T Q = 1//3 C be the op- 
eratin g tem perature. 

In (5.22), with exponentially small error, the main 
contribution to the rj integral comes from the interval 
(0, ju — E + CT ), where E is the bottom of the spec- 
trum and the constant C is several units, Considering 
only band offsets whi ch are linear in the effective mass, 
(|4.39|) , the potential (|J) is 



Ubdd(x) = B[m{x) - moo] + $ sc (x; m; U). (5.23) 
so that 

U v (z) = B[m(z) - moo][l - ij/Bmfz)] + $ sc (x;m; U). 

(5.24) 



We see the condition for neglecting the nonparabolicity 
contribution to the electron density is that (0, fj, — E_ + 
CT a ) is much smaller than the typic al va lues of Bm(z) 
[800-1000meV for Al x Gai- x As, see (jj5|)]. 

After presenting the explicit recipe for computing the 
functional $ sc (z; to; U), let us apply the unitary transfor- 



mation Ux, (L6), associated with the coordinate trans- 
formation (4.3) to the one dimensional Schrodinger equa- 
tion, (4.1) with potential Vs{z). We obtain the BDD- 



type equation (|4.10) with the potential (4.11) and effec- 



tive mass (1121) 

The condition for the two equations to be identical 
is the coincidence of the potentials ( |5.23 ) and U x(x) = 
h 2 Wx{x)/ '2tooo, where Wx(x) is given by ( 4.1l| ). 



Ux{x) = Ubdd(x). 



(5.25) 



This can be rewritten as 



X [x) 1 v[X{x)\ Kqo [ [x>{x)]3 [x , {x)]4 

= -4> sc {x;X';v). (5.26) 



Here. KtQQ IS defined by fl4.42j ) and we introduced the po- 
tentials measured in units Brrioo'- 

v(X) = — ; (j) sc (x;X ;v) = . (5.27) 



Bm^ 



Bm c 



Taking into account ( [4.11| ), ( |4.12| ), ( |4.39| ) and ( |5.16[ ), 

m(x) and Wx (x) are functionals of the coordinate trans- 
formation, X', and t he p otential, Vs, in the original 
Schrodinger equation, ( |4.l| ). We reparameterized the de- 
pendence of & sc (x; to; U) in t erms of these functions, as 
it appears in the second eq. (5-27). 

Setting the self consistency term, 4> sc {x; X'; v), equal to 
zero in (5.26), we r ecove r the differential equation (4.22) 
in the l inear case ( 4. 39] ). The presence of (j> sc (x; X'; v) 
makes ( [5.26 ) a functional equation. 

We will attack this equation in a perturbative man- 
ner, as in section IV E. Indeed, as we have seen there, in 
some interesting cases 1 / Kqo is of the order of the lattice 
spacing. A " quasiclassical" approach to this is to n eglec t 
the term within braces in the right hand side of ( 5.26 ), 
leading to 



X'{x) sa + 



(5.28) 



Although this approximation seems reasonable, we will 
not pursue this further with the same motivation as in 
section [VE. Instead, we will develop again the "small 



potential" approach, using the same motivation as in see- 
Since 4> sc (x; X'; v) is generally smalle r and 



tion 



IV E 



slower varying than i>LY(:e)], the estimate of type ( 4.59] ) 



will be better for 

Before proceeding further, let us write the selfconsis- 
tency term in the right hand side of (5.26) as a function 
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of X(x). As in section [VB, we can t hen reduce the 



order of the derivatives in the equation ( 5.26 ) using the 
fact that it will no longer depends explicitly on x, but 
only through X(x) and its derivatives. Redefining the 
unknown function, the sclfconsistency correction's <f> de- 
pendence on x and the two densities which appear in the 
definition of 4> as functions of X(x), 



S[X(x)]=X'(x); 
4> ac [X(x);S;v] = (j> sc {x; X';v); 

h e i[X(x)] = n e i(x); n d [X(x)] = n d (x). 



and using fl5.8| ) and (5.27), we obtain 

cf> sc (X;S;v) = 



x 



dz 



x 



dt 



X 



dz I dt 

x 



where the we introduced the notation 
27re 2 [h d (z) - n e i(z)\ 



B moc S(z)S(t)e[ moo S 2 (t)Y 



(5.29) 
(5.30) 
(5.31) 



(5.32) 
(5.33) 



It remains to rewrite the density of states v^x^E) as 
a function of X(x). Let us apply the unitary transfor- 



mation Ux to the Schrodinger equation ( |4.l| ) with the 
potential V${z) replaced by 

V Sv {z) = V s {z) + 2m^ [1 - S~ 2 (z)} /h 2 . (5.34) 



The function S(z) was defined in ( 5.29 ). The result of 
the transformation is the BDD equation ( 4.1CQ with the 
potential replaced by 

W Xn (x) = W x (x) + 2 moo r) [1 - S- 2 [X{x)}] /h 2 . (5.35) 



Here Wx is given by ( [4. 11[ ). By ( |5 25 h and the definition 
of S, (5.29), we have that coincides with the 



potential in the equation for the Green function (5.20), 
U v (x). Thus, the Green function G„(x, x', E) is the result 
of the unitary mapping with Ux of the Green function 
of the Schroding er eq uation (4.1 ) wit h potential Vg re- 
placed by Vs n , ( 5.34 ). Finally, ( 5.19 ) can be rewritten 
as 



v v {x,E) = S[X(x)]v 0v (X(x),E). 



(5.36) 



Here, fo?/(z, E) is given by (5.11) in terms of Gori(z, z; E). 
This is the Green functio n of the Schrodinger equation. 
It satisfies the equation (5.12) with m e = moo and the 
potenti al U( z) repl aced by ti 2 Vs r i(z)/2m 00 . Thus, substi- 
tuting (5.36) into (5.22), we obtain the electron density 
as a function of X{x) and a functional of 5* and v. 
We can now rewrite ( |5.26 ) in a form resembling (4.61): 



S" = 3S" 2 /2S + S 2 -l-v + <t> /2, (5.37) 

where we omitted the arguments of the functions. 



The small potential perturbative approach begins with 
introducing the formal small parameter a, setting 



(5.38) 



Then we expand again S int o the serie s (4.60 ) in pow- 
ers of a. We substitute ( 4.6C ) and ( 5.38 ) into ( 5.37 ) and 
expand it into a power series in a. After equating the 
terms with the same power of a and setting a = 1, we 
obtain the "sma ll pot ential" expansion. 

As in section [VE, the zero-th order coincides with a 



uniform structure. The equation S'/ ->(X) = 0, with the 
boundary condition S , ( )(±oo) — * 1, leads to 



5 (0) (X) = 1. 



(5.39) 



The problem with the rest of the series is that although 
the Green function Go can be expanded in terms of the 
potential: 



G = G (0 ) + G (0 )C/G (0 ) + 



(5.40) 



where G^°)(x,y; E) = e ik \ x -y\/2ik is the free particle 
Green function, this expansion assumes that the poten- 
tial is small compared to the kinetic energy. This is cer- 
tainly wrong at the low energies we are interested in. 

We will note instead that the nonparabolicity correc- 
tion to Vs (the second term in ( 5.34 )) is small. To first 
order in t he e xpansion of S in powers of a, using the 
estimate ( |5.28 ), we obtain 



U Sv ~U S + [U s - $ sc ]r)/Bm c 



(5.41) 



Here, Us^ = ft 2 Vsr / /2m oc _and Us — ffVs^m,^ are, re- 
spectively, the potential (5.34) and the reconstructed po- 
tential in the Schrodinger equation measured in energy 
units. We omitted the arguments of the functions. If 
fi — E_ + CT Q is much smaller than Brrioo , then the term 
with square brackets in (5.41) is small compared to the 
first one. 

Neglecting it, the self consistency term coincides in the 
first approximation with the one computed in section V A 
for the Schrodinger dynamics with mass and dielec- 
tric constant e^oo): 

4>{X-l-v) = ^ sc {X;U s )/Bm 00 . (5.42) 

Now, we can solve for the first order correction S\{X): 

1 



Si(X) 



4:K r 



v(0-m\, (5-43) 



with cj) given by (5.42). If the integrand is slowly vary- 
ing on the scale 1/koo, then we recove r the perturbative 
expansion of " quasiclassical" formula ( |5.28 ) . 



2S 1 (X) = v(X)-4,(X;l;v). 



(5.44) 



Final ly, w e find the transformation as at the end of 
section IV E. 
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X(x)=x+ / dySr(y), 
Jo 



(5.45) 



where we inverted the implicit relation ( 4.6§| ) to first or- 
der. The mass (and concentration) profile can now be 
recovered. The result will be reasonable if the electron 
(and dopant) density is not too lar.ge. A more detailed 
account will be published elsewhereES. 



VI. CONCLUDING REMARKS 

We have attained our stated purpose, constructing het- 
erostructure chemical composition and dopant profiles 
such that the electron reflectance of the heterostructure 
is preset, by combining three inverse techniques: 

1. Reconstruction of the phases of the scattering data 
from their absolute values. We did it by Pade dis- 
cretization of the dispersion relations. 

2. Reconstruction of the potential in the one dimen- 
sional Schrodinger equation from the scattering 
data. We used the Kay-Moses-Sabatier approach, 
since the reconstructed scattering data were ratio- 
nal functions of k. 

3. Determination of the chemical concentration and 
dopant ion concentration profiles which yield a 
given self consistent potential for the Schrodinger 
equation or an effective mass profile and a self con- 
sistent potential for the BenDaniel and Duke equa- 
tion. We simply reversed the definition of the self 
consistent potential in the Schrodinger case. In the 
case of BenDaniel and Duke dynamics we sketched 
a perturbative approach to solving the functional 
equation obtained by the variable mass mapping of 
the Schrodinger into the BenDaniel and Duke one - 
the method used for solving the non self consistent 
inverse problem for the BDD equation. 

In section || we used Pade approximation for repre- 
senting the reflectance and reconstructed the phases by 
discrete dispersion relations. The number of parameters 
is determined by two factors: the number of poles of the 
approximants for the transmittance and the reflectance 
(p + q + 2) and the degree of smoothness of the potential. 
The latter is determined by the rate of falloff at infinity 
of the reflection coefficients. If R±(k) ~ k~ q ~ 2 , then the 
q-th derivative of the potential is piecewise continuous. 
The accuracy of representation of the scattering data in- 
creases with increasing p. So does the complexity of the 
calculations. 

The reference solutions are in some sense maximally 
non-symmetric. This can be advantageous in optimizing 
non linear optical response functions, which vanish for 
symmetric potentials. 

We have shown how to construct chemical profiles 
for compositionally graded heterostructures in lattice 



matched systems, like the systems Gai_ x Al x As or 
Ini- x - y Ga x Aly matched on InP, such that in the effec- 
tive mass approximation (with self-consistent potential) 
the conduction electron reflectance (at a preset operating 
temperature) is given by the designer. In the inverse ap- 
proach, the construction of self-consistent potentials can 
be made at little extra cost. Indeed, in the brute- force di- 
rect approach one examines a number of possible config- 
urations and, for each configuration, one has to construct 
iteratively the self-consistent potential and only then one 
can obtain the scattering data for the configuration. 

In the inverse approach for the SE, once the reference 
potential is constructed from the synthesized scattering 
data, one can choose some doping profile and determine 
directly the compositional profile which yields the de- 
sired reflectance. The situation is somewhat more com- 
plicated in the case of the BDD equation. Here, we start 
again from the reference potential. Then, solving a non- 
linear differential equation, we obtained the concentra- 
tion (mass) profile corresponding to a non-selfconsistent 
potential. This is valid if the conduction electron density 
is negligible. 

In the case of non negligible density of conduction 
electrons, a more complicated functional equation must 
be solved. The main difficulties are generated by the 
non-parabolic character of the minibands in the case of 
the BenDaniel and Duke's equation. Perturbative ap- 
proaches were presented in section VB. A hybrid ap- 



proach which looks promising is to search by brute force 
methods near the perturbative solution in the BDD case, 
looks promising. 

We examined only Hartree self consistent potentials. It 
is obvious that the incorporation of exchange-correlation 
corrections which depend only on the local electronic den- 
sity can be done at little extra cost. 

We will mention without details another case in which 
the calculation of the selfconsistent solution is much sim- 
plified. When the electron plasma is thin and hot ev- 



erywhere, e 2 (3 [uei] 1 ^ <C 1 and h 2 (3 [n e i] z/ ^ /m -C 1, we 
may use the (quasi)classical expression for the electron 
density: 

n e i(z; m; U) = n b [m(z)/m 0O ] 3/2 e"™. (6.1) 



2/3 



Throughout this paper we have considered the het- 
erostructure embedded in a homogeneous alloy and did 
not consider the biased case. The methods described 
here can be generalized to the case of different asymp- 
totic compositions at ±oo a«d/or presence of bias. This 
will be presented elsewhereci 

We have examined throughout this paper only con- 
tinuous potentials. If these are digitized using standard 
methods, the low energy dynamics and scattering will 
not be much affected. 

If a more precise, multi-band description is needed, it 
should again be profitable to explore by brute-force stan- 
dard methods the configurations which are close to the 
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ones generated by our one-band approach. Perturbative 
calculations near this point could also be prove useful. 

Finally, let us summarize the advantages of the inverse 
approach over the direct, brute force, one: 

• In practice, the direct method is restricted to opti- 
mizing over small sets of parameters describing the 
structural data. 

• The inverse approach allows the discovery of new 
promising configurations that can be subsequently 
optimized using perturbative and/or traditional 
(brute force) techniques. 

• Obtaining chemical composition and dopant con- 
centration profiles and corresponding to given self- 
consistent potentials for the Schrodinger dynamics 
is rather inexpensive. Although the effort is signif- 
icantly larger in the case of BenDaniel and Duke 
dynamics, good starting configurations can be ob- 
tained perturbatively. 



Let us summarize some properties of the Jost functions 
(JF) which are relevant for the scattering problem: 

a. For real k the JF are continuous in z. / + (z;±fc) 
are a pair of linearly independent solutions of ( |A1| ) . 
The same holds for f~(z; ±fc). 

b. The JF can be continued analytically from the pos- 
itive half-axis (k > 0) to the upper complex half- 
plane Im(fc) > 0. Here, the JF are analytic in k 
with values which are continuous functions of z. 
They have no zeros in z for Im(fc) > 0. 

c. For large complex \k\ in the upper half-plane, the 
JF behave like outgoing waves for all real z: 



f±{z-k) = l + 0{k' 1 )- \k\ 



(A4) 



The above propositions can be readily proved using 
the integral equation of Volterra type which is satisfied 
by the JF: 



APPENDIX A: SCATTERING AND INVERSE 
SCATTERING FOR THE ONE DIMENSIONAL 
SCHRODINGER EQUATION 

For the reader's convenience we will briefly review 
scattering and inverse scattering theory for the one- 
dimensional Schrocfaiges equation. We will use the 
Faddeev-Marchenkcc3~ti2l approach. A detailed account 
may be found in the book by Chadan and SabatierEH. 

Setting the in-plane quasimomentum = brings 



the one-dimensional Schrodinger equation (1.5) to the 
form 



i/>"(z) + [k 2 - V(z)] 1>(z) = 0. 



(Al) 



We will consider only potentials which are piecewise con- 
tinuous and fall sufficiently fast at infinity so that 



dz(l + z 2 )\V(z)\ < oo. 



(A2) 



Then, ( Al) may have only a finite number of bound states 
with negative energy. For E > the spectrum is abso- 
lutely continuous and doubly degenerate. 



1. Jost functions, scattering data and their 
properties 

The scatte ring is best described in terms of the Jost 



solutions of (Al), f±(z;k), which for k = behave like 
outgoing waves near ±oo: 



lim f + (z;k)e 



-ikz 



z — »+oo 



lim f_(z;k)e ikz = 1. 



(A3) 



f ± (z;k)=e ±ikz + 



±oc , sin k(y — z) 

dy ^ L V(y)f±(y,k). 



k 



(A5) 



Now, we have two pairs of linearly independent solu- 
tions of (Al), f±(z; ±k). One of each pair behaves like an 
outg oing/ingoing near the corresponding infinity. Since 
(Al) can have only two linearly independent solutions, 
the outgoing wave JF can be expressed in terms of the 
ingoing wave ones: 

T(k)f + (z; k) = f_(z; -k) + iL (*)/-(*! k); (A6) 
T(fc)/_(z; k) = f + (z; -k) + R+{k)f + (z; k). (A7) 



Instead of seeking the asymptotic behaviors, the trans- 
mission, T(k) and the reflection coefficients to the 
right/left, R±(k), can be expressed in terms of Wron- 
skian determinants of the JF: 



T(k) 



R±(k) 



2ik 



W[U(z;k)J.(z;k)}' 
W[f.(z;±k),f+(z;Tk)] 
W[f+(z;k)J-(z;k)] : 



where we use the notation: 



W[f(z),g(z)} = f'(z)g(z)-f(z)g'(z), 



(A8) 
(A9) 

(A10) 



for the Wronskian of the functions f(z) and g(z). We re- 
mi nd t he reader that the Wronskian of any two solutions 
of ( |Al[ ) does not depend on z and is equal to zero if and 
only if the solutions are linearly dependent. 

For real k the scattering coefficients T(k) and R±(k) 
satisfy the following relations: 
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R+(k)T(-k) + R_{-k)T{k) = 0, 
T{k)T(-k) + R±(k)R±(-k) = 1, 



(All) 
(A12) 



which express the unitarity of the S-matrix. Reality of 
the potential implies also that 



T(k) = [T(-k)]*; R±(k) = [R ± (-k)Y 



(A13) 



Generically, T(0) = and R ± (0) = -1 (if there are 
no "zero energy bound states" - - bounded, but not 
square integrable, solu tions of the Schrodinger equation 
for k = 0). From (A8) and the remark c. above we see 
that T~ 1 (k) can be continued analytically to the upper 
half-plane, lm(k) > 0. Its zeros, if present, are the only 
possible (simple pole) singularities of T[k) in the upper 
half plane. At such a zero, the two JF f±(z; k) are not lin- 
early independent and decay exponentially for z — > ±oo. 
Thus, the poles of T(k) in the upper half plane can occur 
only for k = i\j , where h\j — y/— 2m,, Ki > and Ej are 
the energies of the bound states of (Al). The correspond- 
ing eigenfunctions tpj(z) are real and are normalized by 
J-x dcc|^j (x) | 2 = 1. The two JF are proportional to the 
bound-state wave function: 



f±(z;i\ 3 



= cfM*), 



(A14) 



where Cj are real constants. 



The asymptotic expansion of T{k) near the bound- 



state pole k — Wj is 



T(k) 



A, 



•0(1), 



(A15) 



where C are the constants in (A14). From here we can 



see that the bound-state poles of T(k) must be simple. 
Otherwise, the product of the normalization constants is 
zero. 

Generally speaking, the domains of analiticity of 
R±(k) will be smaller. If the potential V(z) is zero on 
a half-axis and there are no bound states, then the cor- 
responding reflection coefficient is analytic in the upper 
half-plane. Let us check this for the reflection to the 
left coefficient, i?_(fc), and pote ntials which vanish for 
z < 0. The denominator in (A9) is analytic and has no 
zeros since there are no bound states. In the numera- 
tor, f+(z; k) is always analytic and f~(z; —k) = e lkz for 
z < 0. The other case can be dealt with in a similar 
manner. 

From (|A11|-[A1^) we obtain relations between the ana- 
lytic continuations of the scattering data, which are valid 
whenever the arguments of the functions are within the 
domain of analiticity: 

R+(k)T(-k) +R-(-k)T(k) = 0; (A16) 
T(k)T(—k) + R±(k)R±(—k) = 1; (A17) 
[T(k)]*=T(-k*); [R ± (k)}* = R±(-k*). (A18) 

For positive energies, E = h 2 k 2 /2m e , ( A12- A13| ) imply 
that the sum of the transmittance of the heterostructure 
and its reflectance, 



T(E) + 11(E) = \T(k)\ 2 + \R±(k)\ 2 = 1 (A19) 



equals unity. 



2. Inverse scattering. The Marchenko equation. 

If the potential in t he Schrodinger equa tion (|Al[ ) is 
known, then, solving (Al) and using @0) one can 
obtain the scattering data (^2). We want to explore the 
possibility of recovering the pote ntial in the Schrodinger 
equation ( |aT|) from the SD ( ]2~2l) . 

Let 



F ± (z;k)=e* ikz f±(z;k). 



(A20) 



The functions F±(z;k) are analytic in the upper half- 
plane. As a function of k, F±(z; k) — 1 decays at infinity 
no slower than k~ l , (A4). This means that the Fourier 
transformation with respect to k will exist at least in the 
L 2 sense. Wc define the transformation kernels of the 
Schrodinger equation by the Fourier transforms 

1 r +oa 

K + (x;y) = — dfcc- ifc ^)[F + (x;fc)-l], (A21) 

2tt7-oo 

1 f +oa 

K^(x; y) = — dkc ik( - y -^ [F_(x; k) - 1] . (A22) 

27r J-oo 



Closing the integration contour in the upper half-plane, 
we see that 

K + (x,y) = for x > y; K_(x,y) —OiorxKy. 

(A23) 

The inverse Fourier transformations are: 



F+{x;k) = l+ due lk{ - u - x) K+(x;u), (A24) 

J X 

F_(x;k) = l+ [ due ik ( x - u ^K + (x;u). (A25) 



Taking n ow t he Fourier transformations of the Volterra 
equations, (A5), satisfied by the JF yields 



1 [ + °° 
K+{x;y) = - / ds 



V(s) 



+2 



dtV(s-t)K + (s-t;s + t) 



K- (x; y) = - 



V(s) 



+2 



dtV(s-t)K-(s-t;s + t) 



(A26) 



(A27) 
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Finally, taking the limit y — > x ± in (A26-A27) we 
obtain a simple relation between K± and the potential 
V(x): 



Km K+(x;y) = - 

Vlx+O Z 

lim K-(x;y) = - 

yjx — Z 



+ OC 



dsV(s) 
dsV(s). 



(A28) 
(A29) 



Thus, if we can construct K + (or KJ) from_thc scatter- 
ing data, we may use one of the relations ( A28 - A29 ) to 
recover the potential by simple differentiation. 

To make this connection, we will use the analytic prop- 
erties of the scattering data. Applying Cauchy's formula 
to the analytic function F-(x;k) — 1, which decays at 
infinity no slower than fc _1 , we have 



F(x- k) - 1 = lim / dK y ' . 

Sl+o 2tti J_ 00 k — k — id 



(A30) 



Changing k — ► — k in (A30) and using (A7) yields 
1 



F_(x;k) 



1 + lim 

510 27ri /_ 



-°° , 1-T(k)F + (x;k) 
dK — 

K + k + 10 



lim 

540 27ri 



k + k + iS v ; 



In the first integral on the right hand side of (A31) we 
may close the contour in the upper half-plane. If there 
are no bound states, the integrand is analytic in the up- 
per half-plane, and the integral is equal to zero. If bound 
states are present, the first integral in (A31) is equal to 



h 



E 



(Cr) 2 e 2A ^F_(a;;iA J ) 
k + iAj 



(A32) 



where we used (A14) to replace F + (x ; iXA by F_(x; iA 7 ) . 

Taking the Fourier trans form of ( A31 ), using ( |A22| ) 
and the definition, (A25), of K-(x;y), yields the 
Marchenko equation for K_ : 



K-(x;y) + M-(x + y)+ / dsM-(y + s)K^(x; s) = 0. 

J — OC 

(A33) 

Here, the Marchenko kernel, M-(u), is the sum of 
the Fourier transform of reflection to the left coeffi- 
cient i?_(fc) and the contribution from the bound states, 
& 



M_ (u) 



1 

2^ 



/+oo 
dKe-' 1KU R-(n)+J2{Ci 
-CO 



(A34) 



A similar calculation yields the Marchenko equation for 
K+: 



+ 00 



K+(x;y)+M + (x + y)+ / dsM+(y + s)K+(x; s) = 0. 

J X 

(A35) 

The Marc henko kernel M+(u) is given by an expression 
similar to ( A34 ) , with all the — signs changed into + and 
Xj replaced by — Xj. 

The Marchenko equation links directly the SD to the 
transformatio n ker nels K± bypassing the wave functions. 
After solving (A33), the potential V(x) is recovered from 
flA29|) 



V(x) = 2—K_{x;x-0). 
dx 



(A36) 



Let us remar k tha t the first v ariable appears only as a 
parameter in ( A33 ) and-.( A35 ). 

Dcift and TrubowitzO have given necessary and suffi- 
cient conditions for a one-to-one correspondence between 
a set of SD {T(k), R±(k)} whi ch h as no bound states, 
and a potential V(z) sa tisfyi ng (A2) obtained by solving 
Marchenko's equation, (A32). 



1. T(k) and R±{k) satisfy (A11-A13) on the real axis 



2. T(k) is analytic in the upper half-plane and contin- 
uous up to the real axis. 

3. T(k) = 1 + OQkl- 1 ), Im(fc) > 0, and R±(k) = 
0(|fc| _1 ), for real k as |fc| — > oo. 

4- T(k) has no zeros on the real axis, excepting pos- 
sibly a simple one for k = 0. In the latter case, 



l + i?±(0) = 0. 



(A37) 



5. The Marchenko kernels M±(x) are absolutely con- 
tinuous and for any given a > there exists c(a) > 
such that 



dx6[±(x - o)](l + x^)\M' ± (x)\ < c(a) 



(A38) 



The problem becomes a little trickier in the presence 
of bo und states. Then, the Marchenko kernels M±(u), 
( A34 ) , depend now not only on the energy of the bound 
state, which is given by the c orresp onding pole of T(k), 



but also on the constants C , ( A14). The scattering data 



contain information only on the product Cj CJ' , which 



can be recovered from the residue of T(k), (A15). For 
each bound state we can choose one of the parameters 
arbitrarily. The other is fixed by the (analytically 
continued) scattering data. Thus, assuming full knowl- 
edge of the scattering data, in the case when there are 
n > 1 bound states the solution of the inverse problem is 
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not unique. There is a n-parameter family of potentials 
which correspond to the same scattering data. 

From a physical point of view, one cannot recover the 
full information on the bound states in scattering exper- 
iments, which study only the behavior of the solutions 
at large distances, where the relevant information on the 
bound states is exponentially vanishing. 

The numerical solution of the Marchenko equation for 
potential reconstruction is expensive from the computa- 
tional point of view. To find a value for V(x) one has to 
solve ( A33 ) with high enough precision for the subsequent 
numerical differentiation. A lot of useless data is gener- 
ated in th e pro cess, since we need only lim y ^ x K-(x; y) 
for using ( A36 ) (or lim^x K+(x; y)). We will solve the 
Marchenko equation in the manner explained in section 
[II, which is closer to the way we solve the phase recon- 
struction problem in section y. The resource manage^ 
ment compares rather favorably to that of the codesE3'Efl 
which have been written for the direct solution. 



Acknowledgments: We thank Tom Gaylord and 
Elias Glytsis for introducing us to the fascinating sub- 
ject of heterostructure design and for fruitful discus- 
sions. Illuminating discussions with Roger Balian, Pierre 
Sabatier and Giorgio Mantica are gratefully acknowl- 
edged. G.A.M. thanks Alfred Msezane and Carlos Handy 
for support and hospitality at C.T.S.P.S., Clark- Atlanta 
University, where significant parts of the work presented 
here have been done. 



1 G. Bastard, Wave mechanics applied to semiconductor het- 
erostructures; Editions de Physique, Les Ulis, France, 1988. 

2 D.J. BenDaniel and C.B. Duke, Space-charge effects on 
electron tunneling, Phys. Rev., 152, 683 - 692 (1966). 

3 J. Poschel and E. Trubowitz, Inverse spectral theory, Aca- 
demic, New York 1987. 

4 B.N. Zakhariev and A. A. Suzko, Direct and inverse prob- 
lems, Springer, Heidelberg 1990. 

5 B.N. Zakhariev and V.M. Chabanov, New situation in 
quantum mechanics (wonderful potentials from the inverse 
problem), Inverse Problems, 13, R47 - 79 (1997). 

6 V. Milanovic and Z. Ikonic, Equispaced-level hamiltonians 
with variable effective mass following the potential, Phys. 
Rev. B, 54, 1998 - 2003 (1996). 

7 V. Milanovic and Z. Ikonic, Optimization of nonlinear op- 
tical rectification in semiconductor quantum wells using the 
inverse spectral theory, Solid State Commun., 104, 445 - 
450 (1997). 

8 S. Tomic, V. Milanovic and Z. Ikonic, Optimization of 
intersubband resonant second-order susceptibility in asym- 
metric graded Al x Ga±~ x As quantum wells using supersym- 
metric quantum mechanics, Phys. Rev. B, 56, 1033 - 1036 
(1997). 



9 S. Tomic, V. Milanovic and Z. Ikonic, Quantum well shape 
tailoring via inverse spectral theory: optimizing resonant 
second harmonic generation, J. Phys.: Condens. Matter, 
10, 6523 - 6532 (1998). 

K. Chadan and P. C. Sabatier, Inverse Problems in Quan- 
tum Scattering Theory, Springer, New York, 1989. 

11 P. P. Roberts and G. E. Town, Design of microwave fil- 
ters by inverse scattering, IEEE Trans. Microwave Theory 
Tech., 43, 739 - 742 (1995). 

12 D. Bessis, G. Mantica, G.A. Mezincescu and D. Vrinceanu, 
Electron wave filters from inverse scattering theory, Euro- 
phys. Lett., 37, 151 - 155, (1997). 

13 G. Baker jr., Essentials of Pade Approxirnants, Academic 
Press, New York, 1975. 

14 I. Kay, The inverse scattering problem, Research report No. 
EM-74, Inst, of Math. Sciences, Div. of Electromagnetic 
Res., New York University, 1955. 

15 I. Kay and H. Moses, The determination of the scattering 
potential from the spectral measure function III: Calcula- 
tion of the scattering potential from the scattering operator 
for the one-dimensional Schrodinger equation, Nuovo Ci- 
mento, 3, 276 - 304 (1956). 

16 I. Kay and H. Moses, Reflectionless transmission through 
dielectrics and scattering potentials, J. Appl. Phys., 27, 
1503 - 1508 (1956). 

17 I. Kay, The inverse scattering problem when the reflec- 
tion coefficient is a rational function, Commun. Pure Appl. 
Math., 13, 371 - 397 (1960). 

18 P.C. Sabatier, Rational reflection coefficients and inverse 
scattering on the line, Nuovo Cimento, 78B, 235 - 247 
(1983). 

19 R. Balian, D. Bessis and G.A. Mezincescu, Form of kinetic 
energy in effective mass Hamiltonians for hetero structures. 
Phys. Rev. B, 51, 17624 - 17629 (1995). 

20 R. Balian, D. Bessis and G.A. Mezincescu, Mesoscopic de- 
scription of heterojunctions. J. de Physique I, 6, 1377 - 
1389 (1996). 

21 A. Deift and E. Trubowitz, Inverse scattering on the line, 
Commun. Pure. Appl. Math., 32, 121 - 151 (1979). 

22 If TZ(0) < 1, then E — is called a zero-energy bound 



state. The equation (2.1) has a bounded but not square- 
integrable solution. This solution is unstable under small 
perturbations: the zero-energy bound states either disap- 
pears into the continuum or becomes a true bound state. 
In both cases, 1Z(0) = 1. We will not consider the unstable 
case. 

23 W.L. Clinton, Phase determination in X-ray and neutron 
reflectivity using logarithmic dispersion relations, Phys. 
Rev. B, 48, 1 - 5 (1993). 

24 G. Reiss and R. Lipperheide, Inversion and the phase prob- 
lem in specular reflection, Phys. Rev. B, 53, 8157 - 8160 
(1996). 

25 We can regard the eventual double real zeros of V P (E) at 
the transparency energies as con fluen ces of complex con- 
jugate zeros pj ± = p* 2 . Then, in (2.31), the corresponding 
zeros satisfy pj 1 — —p*j 2 - 

26 This smoothness condition ensures that the operator do- 
mains of the initial Schrodinger and the resulting BDD 
type Hamiltonians coincide. In plainer words, if X(x) is 
less smooth, then the equation satisfied by the transformed 



22 



function Uxf will be singular. 

F. Capasso, T. Gaylord and E. Glytsis. Private communi- 
cations. 

G. A. Mezincescu, in preparation. 

D. Bessis, G.A. Mezincescu and P.C. Sabatier. In prepara- 
tion. 

L.D. Faddeev. On the relation between S-matrix and the po- 
tential for the one-dimensional Schrddinger operator, Dokl. 
Akad. Nauk SSSR 121, 63 - 66 (1958); Properties of the S- 
matrix of the one-dimensional Schrddinger equation Trudy 
Mat. Inst. Steklov 73, 314 - 336 (1964). 
V.A. Marchenko, The construction of the potential energy 
from the phases of the scattered waves, Dokl. Akad. Nauk 
SSSR, 104, 695 - 698 (1955). 

V.A. Marchenko, Sturm-Liouville operators and Applica- 
tions, Birkhauser, Basel 1986. 

P.V. Frangos and D.L. Jaggard, A numerical solution to the 
Zakharov-Shabat inverse scattering problem, IEEE Trans. 
Antenn. Propagat., AP-39, 74 - 79 (1990). 
P.V. Frangos and D.L. Jaggard, Analytical and numerical 
solution to the two-potential Zakharov-Shabat inverse scat- 
tering problem, IEEE Trans. Antenn. Propagat., AP-40, 
399 - 404 (1992). 



TABLE I. 12 layer digitized Al c Gai_ c As filter. The Al 
concentrations are c\ — 0.05714, C2 = 2ci, C3 = 4ci; the bulk 
Al concentration is C2 = 11.4%. 



Layer 

# 


Width 
(Atomic layers) 


Width 
(nm) 


Al 

concentration 


1 


6 


1.696 


Cl + C2 + C3 


2 


9 


2.543 





3 


18 


5.088 


C3 


4 


5 


1.413 


C2 


5 


10 


2.827 


Cl + C 2 


6 


10 


2.827 


Cl 


7 


14 


3.957 


Cl + c 2 


8 


14 


3.957 


Cl 


9 


13 


3.675 


Cl + c 2 


10 


13 


3.675 


Cl 


11 


14 


3.957 


Cl + c 2 


12 


11 


3.109 


cl 


FIG. 1. Energy dependence of the transmittance of 12 
layer digitized Al c Gai- c As filter: continuous line — Eq.(4.2); 
dotted line — constant mass approximation Eq.(4.1) with 
mo = m^. Insert: potential energy profile. 




(3)1 
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